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1.  INTRODUCTION 


This  report  is  addressed  to  all  SPICE2  users  involved  in  the  design  of  MOS  (Metal 
Oxide  Semiconductor)  IC’s  (Integrated  circuits).  The  material  contained  herein  serves 
eis  an  addition  to  the  SPICES. G  User's  Guide  [l]. 

The  device  and  model  parameters  documented  refer  to  the  SPICE2.G  releaise  ver¬ 
sions  from  the  University  of  California,  Berkeley  and  obsolete  the  information  con¬ 
tained  in  the  "SPICES  MOS  Modeling  Handbook"  [2]  which  is  vaUd  for  the  SPICES. D 
release  versions. 

The  impact  of  MOS  IC's  in  both  analog  and  digital  appUcations  as  well  as  the 
decreasing  dimensions  of  the  single  transistors  enabled  by  advances  in  processing  have 
made  it  necessary  to  refine  the  models  and  to  provide  more  information  about  each 
device  as  it  appears  on  the  circuit  layout.  Associated  with  each  MOSFET  is  a  drain 
and  source-jimctlon  sidewall  capacitance  (which  has  a  different  voltage  dependence 
than  the  bottom  of  the  diffusion)  and  a  parasitic  series  resistance.  These  are  unique  to 
a  certain  geometry.  At  the  model  level  there  are  effects  which  become  important  as 
the  channel  length  and  width  go  below  IQfxm.  A  thorough  description  of  all  parame¬ 
ters  appearing  on  the  element  (device)  card  and  model  card  is  contained  in  Sec.  2. 

In  SPICE2.G  there  are  three  different  MOS  models  available  to  the  user.  The  Level 
1  model  is  the  simple  Shichman-Hodges  model  implemented  according  to  [3].  This  first 
order  model  has  been  found  necessary  for  checking  out  the  correctness  of  hand  calcu¬ 
lations  when  understanding  or  developing  new  circuits.  The  Level  2  model  is  an  analyti¬ 
cal  one-dimensional  model  which  incorporates  most  of  the  second-order  effects  of 
small-size  devices.  The  Level  3  model  is  a  semi-empirical  model  described  by  a  set  of 
parameters  which  are  defined  by  curve-fitting  rather  than  physical  background.  It  is 
necessary  for  the  circuit  designer  to  know  what  are  the  equations  governing  the 
behaviour  of  the  MOSFET's  and  what  is  the  influence  of  the  various  model  parameters 
even  if  hand  calculations  are  almost  impossible  at  this  complexity.  This  insight  of  the 
M0S2  and  MOSS  model  equations  is  found  in  Sec.  3  and  4. 

The  accuracy  of  the  model  depends  heavily  on  the  values  of  the  input  parameters. 
These  input  parameters  should  be  related  to  the  particular  process  used  at  each 
manufacturing  site.  The  large  number  of  parameters  which  describe  the  model  require 
the  existence  of  a  parameter  extraction  system.  Sec.  5  provides  some  detail  on  how  to 
use  some  simple  preprocessor  programs  for  the  evaluation  of  the  SPICE2  input  param¬ 
eters  from  the  measurements  taken  on  test  devices. 

Sec.  6  gives  an  example  of  a  sample  SPICE2  MOS  IC  ihput  deck  stressing  the  corre¬ 
lation  among  dilTerent  model  parameters  and  the  importance  of  modelling  various 
second-order  effects  for  a  good  agreement  with  measured  pierformance. 
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2.  DESCKIPTION  OF  MDSFOT  PARAMETERS 


2. 1.  Device  Parameters 

According  to  the  new  format  of  the  MOSFET  device  line  [l]  there  are  up  to  8 
pareimeters  describing  each  geometry.  With  the  exception  of  W  and  L  (which  are 
defaulted  to  Im),  all  other  parameters  default  to  zero  if  not  specified. 

W  and  L  are  the  channel  width  and  length  in  meters  as  measured  on  the  layout: 

AD  Eind  AS  are  the  drain  and  source  area,  respectively;  they  scale  the  parameters 
JS  and  CJ  of  the  model  Une.  If  not  specified  they  are  assumed  zero.  When  absolute 
values  for  the  junction  reverse  current  (IS)  or  the  junction  capacitances  (CBD,  CBS)  are 
input  on  the  model  line  the  areas  can  be  omitted.  AD  and  AS  are  measured  as  if  on  a 
layout  and  thus  model  the  bottom  of  the  junction, 

PD  and  PS  are  the  perimeters  of  the  drain  and  soiu-ce  respectively;  they  multiply 
the  junction  sidewall  capacitance  CJSW  specified  in  F/m  on  the  model  line. 

NRD  and  NRS  are  the  number  of  squares  of  the  parasitic  series  resistance  of  the 
drain  and  source  diffusions  as  estimated  from  the  layout.  NRD  and  NRS  multiply  the 
sheet  resistance  RSH  given  on  the  model  card.  An  alternate  way  to  specify  series  resis¬ 
tances  for  the  drain  and  source  is  to  enter  the  total  values  RD  and  RS  on  the  model 
card;  these  values  will  then  be  added  to  all  devices  which  invoke  that  particular  model. 
When  RD  and  RS  appear  in  the  definition  of  a  model  then  those  devices  referencing  this 
model  need  no  NRD  or  NRS. 

As  is  seen  from  the  above  there  is  a  close  correlation  between  the  parameters 
which  appear  on  the  device  card  and  the  parameters  input  on  the  model  card.  Since  a 
hiereirchical  order  of  the  model  parameters  is  given  in  the  following  paragraph, only  a 
single  example  is  mentioned.  The  presence  of  NSUB  on  the  model  line  automatically 
calls  for  the  calculation  of  CJ.  A  way  to  strip  the  model  of  cuiy  bottom  junction  capaci¬ 
tance  is  to  delete  AD  and/or  AS  from  the  device  card. 


2.2.  Model  Parameters 

This  section  rather  than  reproducing  the  information  contained  in  the  User’s 
Guide  groups  the  37  parameters  describing  the  MOSFET  model  in  a  logical  meinner. 

LEVEL  indicates  the  program  which  model  is  desired; 

LEVEL=1  invokes  the  Shichman-Hodges  model  (default); 

LEVEL=2  invokes  the  "M0S2"  model  to  be  described  in  the  following  section; 

LEVEL=3  invokes  the  semi-empirical  model  described  in  Sec.  4. 

A  first  classification  is  to  di'dde  the  input  data  into  "Electrical"  (derived)  and  "Pro¬ 
cessing"  (primary)  parameters.  According  to  this  view  point,  VTO  and  GAMMA  are 
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electrical  and  NSUB  and  TOX  are  processing  peu’ameters.  For  all  models  both  kinds  of 
parameters  can  be  entered  with  the  general  convention  that  the  electrical  data  will 
always  override  the  value  computed  from  processing  data  if  edso  specified.  Thus  if  VTO. 
NSUB  and  TOX  are  input  the  threshold  voltage  wUl  assume  the  value  entered  as  VTO, 
while,  e.g.,  GAMMA  will  be  computed  from  NSUB  and  TOX. 

The  following  5  data  are  the  Electrical  (derived)  parameters; 

VTO  is  the  Extrapolated  Zero-Bias  Threshold  Voltage  of  a  long  and  wide  channel 
device.  This  voltage  is  derived  at  the  onset  of  strong  inversion  and  marks  the  point 
where  the  device  starts  conducting  if  the  weak  inversion  current  is  neglected  (see  the 
description  of  the  parameter  NFS). 

If  VTO  is  not  an  input  parameter  and  NSUB  and  TOX  are  specified,  SPICE2  will 
evaluate  the  threshold  from  the  equation 


Vto  -  VpB  +  Ztpf  + 


Cox 


(2.1) 


where 


Yfb  -  <pys  - 

is  the  flatband  voltage  [4],  [5]  and  is  the  meted  (poly-silicon)  semiconductor  work 
function  difference.  VTO  will  be  printed  out  in  the  summary  of  the  model  parameters 
at  the  beginning  of  the  SPICE2  output.  The  value  VTH  which  is  printed  in  the  operating 
point  information  represents  the  threshold  for  given  terminal  voltages  and  for  the  par¬ 
ticular  device  including  size  dependence,  (see  Sec.  3.1.). 

KP  is  the  Intrinsec  Transconductance  Parameter.  If  it  is  not  specified  and  UO  (sur¬ 
face  mobility)  and  TOX  are  entered,  it  is  computed  by  the  program  according  to  the 
equation 


KP-UOy-Cax 


(2.3) 


GAMMA  is  the  Bulk  Threshold  Peirameter  and  represents  the  proportionality  factor 
relating  the  change  in  threshold  voltage  to  backgate  bias.  It  is  encountered  in  the 
expression  of  \T0,  the  dreiin-source  current  lss<  and  saturation  voltage.  When  it  is  not 
input  and  if  the  necessary  processing  parameters  are  given,  it  is  derived  from 


GAMMA  = 


V2g  ssiMSUB 

Cox 


(2.4) 


if  the  necessary  processing  parameters  are  given. 
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PHI  is  the  Surface  Potential  at  strong  inversion  2^5 When  not  input  it  is  com¬ 
puted  from 


PHI  =2 


kT  ^^NSUB 
q  ni 


(2.5) 


LAMBDA  is  the  Channel-Length  Modulation  Parameter  and  is  equivalent  to  the 
inverse  of  the  Early  voltage  for  the  bipolar  transistor  and  is  a  measure  of  the  output 
conductance  in  saturation.  By  specifying  this  parameter  the  MOSFET  will  have  a  finite 
but  constant  output  conducteuice  in  saturation. 

If  LAMBDA  is  not  input  the  LEVEL=1  model  will  assume  a  zero  output  conductance 
while  the  LEVEL=2  model  will  compute  a  fi^nite  and  voltage-dependent  output  conduc¬ 
tance  defined  by 


LAMBDA  = 


LL 

LxVds 


(2.6) 


where  L  is  the  channel  length  input  on  the  device  line  less  the  lateral  diffusions  of  the 
drain  and  source,  LD,  and  AX  is  a  function  of  Vos  and  Vqsat  is  shown  later. 

The  above  set  of  electrical  parameters  is  adequate  for  the  LEVEL=1  model. 

In  the  following  the  parameters  which  characterize  the  gate  and  channel  of  the 
MOSFET  are  summarized.  These  parameters  are  further  classified  as  oxide  charac¬ 
teristics.  charge  concentrations  and  mobility  parameters. 

The  characteristics  of  the  thin  oxide  are  listed  below. 

TOX  is  the  Oxide  Thickness  and  comes  into  the  calculation  of  the  conduction  fac¬ 
tor,  backgate  bias  effect  and  gate-channel  voltage-dependent  capacitances. 

If  not  specified,  in  the  1jEVEL=1  model  TOX  will  be  assumed  infinite,  i.e.,  VTO,  KP,  or 
GAMMA  will  be  defaulted  rather  than  computed  if  these  last  quantities  are  not  entered 
directly.  The  gate  capacitances,  Cos,  Cqo  and  Cqb  are  assumed  to  be  the  overlap  capa¬ 
citances,  i.e.,  they  are  constant  with  voltage.  In  the  liEVEL=2  and  =3  model,  TOX 
defaults  to  lOOOA  implying  that  the  gate  capacitances  are  always  voltage-dependent. 

CGSO  is  the  Gate-Source  Overlap  Capacitance  per  meter  channel  width  and 
represents  the  capacitance  due  to  the  lateral  diffusion  of  the  source  in  a  silicon  gate 
MOSFET;  the  input  value  is  multiplied  by  the  channel  width.  This  is  the  constant  part  of 
the  total  Cos  capacitance  which  has  also  a  contribution  due  to  the  thin  oxide  capaci¬ 
tance  (the  voltage-dependent  charges  associated  wth  the  gate  and  channel). 

CGDO  is  the  Gate-Drain  Overlap  Capacitance  per  meter  channel  width.  The  same 
considerations  as  above  apply. 

CGBO  is  the  Gate-Bulk  Overlap  Capacitance  per  meter  channel  length  and  is,  as 
opposed  to  the  above  capacitances,  a  field  oxide  capacitance.lt  is  the  result  of  the 


-4- 


requirement  that  the  gate  has  to  extend  beyond  the  channel  by  a  certain  amount 
according  to  the  specific  design  rules.  The  total  Cgp  overlap  capacitance  results  by 
multipljhng  the  above  value  by  the  channel  length. 

There  are  five  parameters  which  describe  or  relate  to  charge  concentrations. 

NSUB  is  the  Substrate  Doping  and  is  used  in  the  derivation  of  most  of  the  electri¬ 
cal  parameters.  The  effect  of  implants  can  be  added  by  specif3dng  different  VTO's  and 
maintaining  the  same  NSUB  in  the  different  models  for  the  calculations,  say  of  the 
junction  capacitances.  NSUB  can  be  specified  for  a  LEVEL=1  model  as  well  and  it 
assumes  the  same  role. 

NSS  is  the  Effective  Surface  Charge  Density  and  is  used  for  the  evaluation  of  VTO. 
In  most  cases  it  is  more  accurate  to  specify  \TO  as  obtained  from  a  parameter  extrac¬ 
tion  system  rather  than  NSS. 

NFS  is  the  Effective  Fast  Surface  State  Density  and  serves  both  as  a  flag  for  the 
evaluation  of  subthreshold  currents  and  as  a  quantity  controllir-g  the  amount  of 
subthreshold  current  flowing  in  the  device.  If  not  specified,  the  MOSFET  conducts  only 
in  strong  inversion.  This  parameter  is  meaningless  for  a  LEVEL=1  model. 

NEFF  is  the  Total  Channel  Charge  (Fixed  and  Mobile)  Coefficient.  It  is  used  as  a 
multiplicative  factor  of  NSUB  in  order  to  get  the  proper  value  of  the  output  conduc¬ 
tance  in  saturation.  Physically  it  accounts  for  the  fact  that  in  the  channel  there  is 
both  mobile  and  fixed  charge  and  the  total  amount  has  to  be  considered  for  computing 
the  output  conductance.  NEFF  has  meaning  only  when  the  velocity  saturation  model  is 
used  (VMAX  is  specified)  in  the  LEVEL=2  model. 

TPG  is  the  Tjq>e  of  Gate  indicates  whether  the  simulated  device  has  a  metal  or 
polysilicon  gate  (and  which  polarity  type  relative  to  the  substrate).  The  NGATE  param¬ 
eter  of  the  older  SPICES  (D  and  E)  versions  has  been  deleted  since  the  poly  gate  is  usu¬ 
ally  heavily  doped  and  the  Fermi  level  is  located  very  close  to  one  of  the  band  edges. 

A  last  group  of  parameters  characterizing  the  gate  and  channel  refer  to  the  pro¬ 
perties  of  carriers  in  the  conductive  channel.  The  first  foui'  describe  the  surface  mobil¬ 
ity  and  its  variation  with  the  electric  field. 

UO  is  the  Surface  Mobility  at  low  gate  voltages.  It  is  the  only  mobility  parameter 
which  has  meaning  for  a  LEVEL=1  or  =3  mcdel. 

UCRIT  is  the  Critical  Field  for  Mobility  Degradation  and  is  the  limit  at  which  the 
surface  mobility  UO  starts  decreasing  according  to  the  empirical  relation  given  below. 

UEXP  is  the  Critical  Field  Exponent  for  the  empirical  formula  which  characterizes 
the  degradation  of  the  surface  mobility. 

UTRA  is  the  Transverse  Field  Coefficient  for  the  empirical  mobility  degradation  for¬ 
mula. 
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The  above  four  parameters  enter  in  the  following  equation  of  the  surface  mobility 
degradation  [6] 


IJ.S  =  I70x 


_ UCHITxssi _ 

Coxi^Gs  —  Vm  —  UTEAxVos) 


VEXP 


(2.7) 


The  above  parameters  are  determined  best  by  curve  fitting  as  is  shown  in  Sec.  5. 
In  most  cases  UTRA  can  be  chosen  between  0.0  and  0.5. 

VMAX  is  the  Maximum  Drift  Velocity  of  Carriers  and  affects  the  drain-source 
saturation  voltage  of  the  MOSFET.  It  has  been  proven  that  for  devices  with  channel 
lengths  shorter  than  10ft  the  drain  current  saturates  before  the  channel  pinches  off 
because  the  electrons  (holes)  reach  the  scattering  limited  velocity  in  the  channel.  The 
model  proposed  by  Baum  and  Beneking  [7]  and  first  used  in  the  SHOMOS  model  [8]  is 
applied  in  order  to  get  Vq^at  from  the  following  equation 


VMAX  = 


IpSAT 

^’^Qchan 


(2.8) 


where  Idsat  Qchan  functions  of  Vjjsat-  More  details  about  this  are  given  in  the 
following  section. 

VMAX  also  acts  as  a  flag  in  the  l.EIVEL=2  model.  When  specified  the  scattering  lim¬ 
ited  velocity  of  carriers  determines  Vosat\  otherwise  the  pinch-off  approach  gives  Vdsat- 
This  parameter  is  meaningless  for  the  LE\'EL=1  model.  When  specified  in  LEVEL=2,  the 
channel  shortening  effect  uses  the  NEFF  parameter  as  introduced  above  (see  also  Sec. 
S). 

A  last  category  of  parameters  characterize  the  drain  and  source  junctions  of  the 
MOSFET.  There  is  some  overlap  among  some  of  these  parameters,  e.g..  the  reverse 
current  of  the  j’jnction  can  be  input  either  as  IS  (in  A)  or  as  JS  (in  A/m,®).  Whereas  the 
first  is  an  absolute  value,  the  second  is  multiplied  by  AD  and  AS  to  give  the  reverse 
current  of  the  drain  and  source  junctions  respectively.  This  methodology  has  been 
chosen  since  AD  and  AS  default  to  zero.  The  same  procedure  applies  also  to  the  zero- 
bias  junction  capacitances  CBD  and  CBS  on  one  hand,  and  CJ  on  the  other  as  is  seen  in 
the  following  description. 

RD  andRS  are  the  Drain  and  So'irce  Ohmic  Resistances  and  are  values  which  apply 
to  all  devices  invoking  the  model  in  which  they  are  specified.  (This  was  the  only  way  for 
the  SPICE2.D  and  .E  versions  to  include  series  resistances). 

RSH  is  the  Sheet  Resistance  of  the  Drain  and  Source  Diffusions  and  is  multiplied  by 
the  number  of  squares  NRD  and  NRS  which  are  entered  on  each  device  card  to  give  the 
total  drain  and  source  ohmic  resistcuice  respectively,  for  each  device. 
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When  both  the  total  value  (RD  or  RS)  and  the  RSH  are  input  the  total  value  over¬ 
rides. 

CBD  and  CBS  are  Zero-Bias  B-D  and  B-S  junction  capacitances  and  are  total  values 
for  all  devices  which  reference  the  model  where  they  are  specified. 

CJ  is  the  Zero-Bias  B-D  and  B-S  junction  capacitance  per  of  junction  area  and  is 
multiplied  by  the  areas  AD  and  AS  entered  on  the  device  card  to  give  total  values  for 
the  B-D  and  B-S  junction  capacitances,  respectively,  specific  to  each  single  device.  CJ 
is  computed  automatically  by  the  program  if  NSUB  is  specified  and  a  step  junction  is 
assumed. 

When  both  the  total  value  (CBD  or  CBS)  and  CJ  (br  NSUB)  are  input,  the  total  value 
overrides. 

MJ  is  the  Grading  Coefficient  of  the  B-D  and  B-S  junctions  eind  applies  to  both  CBD 
euid  CBS.  or  CJ. 

The  above  capacitances  are  well  suited  for  modeling  the  bottom  of  the  junctions. 
It  has  been  found  experimentally  that  an  important  part  of  the  junction  capacitance  is 
contributed  by  the  sidewall  of  the  diffusion  vdiich  extrapolates  to  a  different  zero-bias 
value  and  grading  coefficient  than  those  pertaining  to  the  bottom  of  the  diffusion.  The 
followixig  two  parameters  account  for  this  capacitance. 

CJSW  is  the  Zero-Bias  Junction  Sidewall  Capacitance  per  meter  of  Drain  and  Source 
Perimeter.  It  is  multiplied  by  PD  and  PS  entered  on  the  device  line  to  give  the  drain 
and  source  junction  sidewall  capacitances,  respectively. 

MJSW  is  the  Grading  Coefficient  of  the  Junction  SidewaE  . 

All  the  above  capacitances  and  grading  coefficients  are  determined  best  from 
curve  fitting  of  measurements  taken  on  specially  designed  test  chips  as  is  outlined  in 
Sec.  5. 

The  total  junction  capacitance  cam  be  expressed  as  follows: 

For  reverse  bias,  Vgs  <  FCxPB  (see  below  for  explanation  of  FC  and  PB) 


CBSf(j(r  ~  CJ 


AS 


^BS 

Tsr 

[  ^ 

+  CJSW 


PS 


1  - 


PB 


(2.9) 


For  forward  bias,  Vgs  >  FC^-PB 

AS 


CBSj’fff  —  CJ 


(1 


1  -  FC{1+MJ)  + 


(2.10) 


-7- 


CJSW 


PS 

(1  — 


1  -  FC{l-i-MJSW)  + 


Vbs 

PB 


^MJSW\ 


The  above  equations  insure  a  smooth  transition  from  the  reverse  bias  region  to  the 
forward  bias  region  and  prevent  the  capacitance  from  increasing  to  infinity 

IS  is  the  Reverse  Current  of  the  Drain  or  Source  Junctions  and  is  a  total  value,  the 
same  for  all  the  drain  and  source  junctions  of  the  devices  which  reference  the  model 
containing  IS, 

JS  is  the  Reverse  Current  Density  of  the  Drain  or  Source  Junction  per  of  junc¬ 
tion  area.  The  total  value  IS  results  from  a  multiplication  by  AD  emd  AS  and  is  specific 
to  each  junction  of  each  transistor. 

If  neither  of  the  above  currents  are  specified,  the  program  will  default  IS  to  l.OE- 
14  A.  If  both  IS  and  JS  are  specified,  IS  overrides  JS. 

PB  is  the  Bulk  Junction  Potential  is  used  in  the  junction  capacitance  formulation. 

FC  is  the  Forward  Bias  Nonideal  Junction  Capacitance  Coefficient  and  is  used  as 
shown  above  in  matching  the  transition  point  and  voltage  characteristic  of  the  junction 
capacitance  when  bias  changes  from  reverse  to  forward.  The  default  value  of  0.5 
proves  satisfactory  most  of  the  time. 

XJ  is  the  Metallurgical  Junction  Depth, 

LD  is  the  Lateral  Diffusion  Coefficient  in  meters.  The  effective  channel  length  is 


L  —  Ly  —  Zy-LD 


(2.11) 


where  Ly  is  the  length  on  the  circuit  mask  input  on  the  device  line.  In  all  subsequent 
equations  L  represents  the  effective  channel  length. 

LD  is  used  for  the  computcition  of  L  only,  whereas  XJ  is  both  a  flag  and  parameter 
used  in  the  evaluation  of  the  short-channel  effect. 

The  only  parameters  not  mentioned  so  far  are 

KI^  is  the  Flicker  Noise  Coefficient 

AF  is  the  Flicker  Noise  Exponent.These  two  parameters  are  used  in  the  small- 
signal  AC  noise  analysis  to  determine  the  equivalent  noise  ciorrent  generator  connected 
between  dreiin  and  source,  the  value  of  which  is 


KFxI;^ 
f  Cox  W  L 


In  the  future  an  improved  noise  model  for  the  MOSFET  is  needed. 


(2.12) 
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Listed  below  are  four  model  parameters  which  are  specific  to  the  LEVEL=3  model 
only.  For  this  model,  the  parameters  LAMBDA,  UCRIT,  UEXP,  UTRA,  and  NEFF  are  mean¬ 
ingless. 

ETA  is  the  Static  Feedback  Effect  Parameter.  It  is  the  proportionality  factor  which 
multiplies  the  value  of  1^5  in  the  Vth  equation  implemented  in  M0S3  to  model  the 
static  drain  to  gate  feedback  effect. 

DELTA  is  the  Channel  Width  Factor.  This  is  an  empirical  factor  which  adjusts  the 
depleted  bulk  charge  underneath  the  gate  to  be  in  agreement  with  the  observed  thres¬ 
hold  voltage  increase  in  narrow-channel  devices. 

THETA  is  the  Empirical  Mobility  Modulation  Parameter  and  is  used  to  compute 
mobility  degradation  as  a  function  of  Vqs- 

KAPPA  is  the  Field  Correlation  P’actor.  It  is  used  to  adjust  the  electric  field  across 
the  surface  depletion  region  of  the  device  in  satiaration. 
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3.  IlCTJEmNTATION  OF  SHAUrSIZE  DEVICE  SECONDORDER  EFFEECTS  IN  THE 

MOS2  MODEL 

This  section  presents  the  main  features  of  the  LEVEL=2  (M0S2)  model  and  the 
equations  which  describe  the  various  first-  and  second-order  effects. 

The  following  effects  are  modeled  in  M0S2; 

-  Backgate  bias  and  small-size  (short-  and  narrow-channel)  effect  on  threshold  voltage 
considering  also  the  static  feedback  from  drain  to  gate; 

-  Saturation  due  to  the  scattering  limited  drift  velocity  of  carriers  and  finite  voltage- 
dependent  output  conduct£mce; 

-  Surface  field  dependent  mobility: 

-Weak  inversion  conduction; 

-  Charge  controlled  mociel  of  regenerative  effects; 

-Variation  of  all  quantities  with  temperature. 

In  the  rest  of  this  section  the  pertinent  equations  are  listed  with  emphasis  on  the 
influence  of  the  different  input  parameters. 

3.1.  Threshold  Voltage 

From  the  zero-bias  threshold  voltage  VTO  given  in  Sec.  2.  which  applies  for  a  MOS- 
FET  with  large  dimensions  (W,  L  >  20/2m)  and  with  the  source  and  bulk  connected 
together  SPICE2.G  readjusts  this  value  to  Yfu  at  each  operating  point  and  for  each 
transistor  according  to  its  size. 

Physically  the  following  effects  take  place; 

-  An  increase  of  the  bulk-to-source  voltage  increases  the  depletion  charge  which  causes 
VfH  to  go  up. 

-  In  a  short-channel  device  part  of  the  depletion  charge  in  the  bulk  terminates  the  elec¬ 
tric  field  of  the  drain  and  source  junctions.  The  electric  field  from  gate  to  bulk 
depletes  less  charge  and  thus  Vth  is  lowered.  This  can  be  viewed  by  the  trapezoid 
approach  introduced  by  Yau  [9]  and  which  is  used  by  SP1CE2  (see  Fig.  3.1.). 

-  The  edge  effects  in  a  narrow  channel  cause  the  depletion  charge  to  extend  beyond  the 
width  of  the  channel.  The  gate-to-bulk  field  has  to  be  increased  to  balance  this  charge; 
thus,  V-pH  increases  (see  Fig,  3.2,). 

-  The  amount  of  charge  underneath  the  gate  depleted  by  the  drain  junction  field  rather 
than  the  gate-to-bulk  field  increases  with  ^DS  thus  lowering  This  effect  is  approxi¬ 
mated  geometrically  in  M0S2  as  shown  in  Fig.  3.3.  and  is  referred  to  as  the  static 
drain-to-gate  feedback  on  threshold. 

An  equation  for  VpH  which  sums  up  the  above  features,  can  be  expressed  as 
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(X+XJ)^  =  (XJ-Wq)^  - 
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FIG.  3.1, 
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FIG.  3.2. 
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Vf}j  -  Vgju  7sy/2(pp  -  Vgs 


(3.1) 


where 

^BiN  =  ^/  +  DELTA-^^^^{2(Pp-Vbs)  (3.2) 

is  the  corrected  built-in  voltage  for  narrow  channel,  and 


(3.3) 


75  is  the  corrected  GAMMA  for  short-chauinel  with  static  drain-to-gate  feedback  (see 
Fig.  3.3.) 


75  =  GAMMA  {l-OLs-oii)) 


(3.4) 


where  and  are  the  correction  factors  for  the  depletion  charge  at  the  drain  and 
source,  respectively, 


“5 


1  XJ 

2  L 


H 

XJ 


1 


(3,5) 


O-D 


1  XJ 

2  L 


1 


(3.6) 


The  depletion  widths  Ws  and  Wj)  are 


-^D  '^^9f~^bs 


(3,7) 


Wj}  —  Xo'\/Zff  —  Vqs  4  Vds  (3.8) 

where 


3a  SI 
NS  US 


(3.9) 


The  implementation  of  the  geometry  in  Fig,  3.3.  represents  an  idealization  of  the 
actual  depletion  layer  shape.  Preliminary  results  show  that  it  offers  a  good  approxima¬ 
tion  at  a  minimum  computational  effort. 

It  should  be  noticed  that  the  absence  of  XJ  from  the  input  file  bypasses  the  correc¬ 
tion  for  short  channel  effect. 
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3.3.  Dram-Source  Current 

In  order  to  correct  for  small  dimensions,  the  implemented  drain-source  current 
equation  is  slightly  modified  relative  to  previous  formulations  [4],  [5]  which  include  the 
effect  of  the  bulk  charge. 

The  current  in  strong  inversion  ( Vcs  >  V^f)  is  given  by 


hs  - 


Vqs  -  Vgjff  - 


V^DS 


VdS  -  ^s\'ZfF^Vos-VBS?'^  - 


(3,10) 


where 


7]=!^  DELTA 


ACoxW  ^ 


^--J^sCqX 

It  can  be  seen  that  the  nairow  channel  effect  is  included  through  Vbib  given  in  Eq. 
(3.3)  and  r].  The  short  channel  uses  the  corrected  GAMMA.  75,  ns  represents  the 
degraded  surface  mobility.  For  large  "Vf  and  L  the  above  equations  revert  to  the  well 
known  forms  [4].  [5]. 


3.3.  (^ration  in  the  Saturation  Region 


3.3.1.  Saturation  due  to  Pinch-off;  Channel  length  Modulation 

As  mentioned  earlier,  if  VMAX  is  not  input  the  program  computes  the  saturation 
voltage  Vbsat  assuming  the  channel  pinches  off  at  the  drain  limit.  With  the  corrections 
for  small-size  effects.  1^547’  is 


''DSAT  - 


{Vcs-Vblw)  .  1 

7s 

77  3 

.  'n  . 

>2 


1+4 


V 

75  J 


Vcs  -  ^BIN 


i-2cpB~^BS 


1/2 


(3.11) 


The  finite  output  conductance  observed  in  saturation  is  due  to  the  spread  of  the 
pinch-off  region  into  the  channel  thus  reducing  the  effective  channel  length. 
L,ff=L—LL,  and  increasing  the  geometrical  ratio  of  the  conduction  factor. 


W  _ W _ 

L  -M.  L{1-  LAMBDAx  V^s) 


(3.13) 
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LAMBDA  =  ,  „ 
LxVjos 


If  LAMBDA  is  not  an  input  parameter,  its  value  is  computed  according  to  Eq.  (3.12) 
where 


LL=Xo 


yOS-Vl^AT 


V 


1+ 


yiS~yDSAT 


1/2 


(3.13) 


This  formulation  does  not  include  the  effect  of  the  field  between  gate  and  drain  and 
gate  and  pinch-off  point  respectively  [10],  but  insures  the  continuity  of  the  current  and 
its  first  derivative  at  the  transition  from  the  triode-  into  the  saturation  region  (impor¬ 
tant  in  an  iterative  solution  algorithm  of  nonlinear  simultaneous  equations).  Another 
drawback  of  the  above  formulation  of  the  chamiel  modulation  is  that  it  overestimates 
the  output  conductance  in  saturation. 

One  problem  with  short  channel  devices  is  that  at  high  Vj}s  punch-through  occurs. 
This  effect  is  not  modeled  in  SPICE2.  In  order  to  prevent  the  channel  length  from  going 
negative,  Eq.  (3.12)  is  used  only  when  the  effective  channel  length  is  larger  than 
the  zero-bias  depletion  layer  width  Wb,  {WB=XD\rFB).  When  L^jj  gets  smaller  than 
Wb,  Eq.  (3.12)  is  replaced  by 


A.//  = 


^'b 


1+- 


Wb 


(3,14) 


where 

=  L-Wb  (3.15) 

Although  this  equation  does  not  model  punch-through  it  prevents  numerical  non- 
convergence  problems  that  have  been  observed. 


3.3.2.  Saturation  due  to  Scattering  Limited  Velocity,  Channel  Length  Modulation 

In  short-channel  MOSFET's  the  saturation  of  the  current  occurs  because  the 
charge  carriers  reach  their  maximum  scattering  limited  velocity  before  pinch-off.  In 
other  words,  a  lov/er  current  is  expected  from  a  short  channel  device  as  compared  to  a 
long  channel  device  of  identical  geometric  ratios,  processing  and  biasing  conditions. 

Vbsat  is  obtained  as  the  solution  of  Eq.  (2.8)  after  substituting  the  expressions  for 
^DSAT  Qchan-  This  leads  to 
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f^s 


VMAX  = 


■  Vos’  -  VbIN  -  ^  VdSAT  -  l^s  (  [vb5.4r +2i;Sf  - 1^5  j  ®  V55  j  ® ) 


J^8ff{yGS-yBIN^VDSAT-7s\yDSAT  +  ^‘pF-VBS^  ) 


.(3.16) 


A  problem  arises  because  Qcban  is  defined  at  the  point  where  the  carriers  reach  their 
maximum  velocity.  This  point  moves  towards  the  source  as  V^s  increeises.  The  shorten¬ 
ing  of  the  channel  is  implemented  in  M0S2  as  proposed  by  Baum  [7], 


L,ff 


=  L-Xo-\/ 


Xd  VMAX 


2^5 


+  {VDs-VosAr)+ 


X§  VMAX 
Zfis 


(3.17) 


In  this  theory  on  a  Vcs-=constant  curve,  the  saturation  voltage  varies  (decreases)  as 
Vps  increases.  From  a  numerical  point  of  view  this  requires  the  simultaneous  solution 
of  two  nonlinear  equations,  i.e.,  iterative  computation  which  increases  the  computer 
analysis  time.  There  is  a  proposal  [11]  which  limits  the  solucion  to  two  iterations.  How¬ 
ever  the  method  adopted  in  SPICES. G  precludes  any  nested  iterative  process  based  on 
the  following  assumptions: 

First,  define  Vj)sat  according  to  Eq.  (3.16)  with  Lgff=L  and  use  thj.s  value  for  any  Vos  on 
a  FciS=  constant  characteristic.  This  is  based  on  the  fact  that  Eq.  (3.17)  does  not  pro¬ 
vide  an  accurate  description  of  the  output  conductance  in  saturation  and  NEFF  has  to 
be  used  as  an  empirical  factor.  In  the  expression  of  Xo,  Eq.  (3.9),  NSUB  is  replaced  by 
NEFF*NSUB. 

Second,  good  agreement  of  the  simulated  and  measured  1-V  characteristics  can  be 
obtained  by  adjusting  NEFF. 

A  closed-form  solution  can  be  obtained  as  follows.  Eq.  (3.16)  can  be  simplified  to 
obtain 


V 


^{x^-Vz)-^s{^^-Vo^) 
Fi— 75Z  —x^ 


(3.18) 


where 


-  >/  Vdsat'^'^Vf~Vbs 
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VMAX^L.ff 

v  = - 

t^s 

Eq.  (3.16)  can  be  written  as  a  quartic  equation 
x*+ax^-i-bx^+cx  +d  =0 
with 


(3.19) 


fe=-2(Ki+u) 

c=-2ys'^ 


<i  =  -V2^-^sy2^+^V,V^^2vV, 

This  equation  can  be  solved  by  Ferrari’ s  method.  Among  the  2  or  4  real  roots,  the 
smallest  positive  one  is  the  valid  solution:  the  choice  is  based  on  the  fact  that  the  equa¬ 
tion  represents  the  difference 


fDSAT~y^AX'!<WxOcHAN{i^)  (3.20) 

which  decreases  as  Vus^T  increases. 

In  some  cases  the  equation  might  not  have  a  real  solution.  In  this  case  the  pro¬ 
gram  computes  V^sat  based  on  the  pinch-off  approximation  and  prints  out  a  message. 


3.4.  Weak  Inversion  Conduction 

A  MOSFET  is  not  an  ideal  svhtch  which  starts  conduction  abruptly;  there  is  current 
flowing  in  the  device  below  the  threshold  voltage  which  marks  ocily  the  onset  of  strong 
inve.rsion.  This  current  is  known  as  the  weak  inversion  or  subthreshold  current.  The 
simvilation  of  this  behaviour  is  important  for  present  day  MOS  circuits  which  eire  often 
de.signed  to  operate  in  the  weak  inversion  region.  The  model  implemented  in  SPICE2.G 
uses  as  a  starting  point  the  analysis  performed  by  Sv/anson  and  Meindl  for  CMOS  dev¬ 
ices  [12]. 

As  shown  in  Fig.  3.4,  a  new  threshold  voltage  Vqv  is  defined,  above  %/.  which 
marks  the  transition  from  the  weexk  inversion  to  the  strong  inversion  characteristic; 


(3.21) 


nkT 

1 


where 


n  =  1  + 


(^FS 

Cox 


Cp 

Cox 


(3.22) 


Cps  =7  ^NFS 


(3.23) 


_  ^Qb  _ 
~  dVBs  ~ 


d-\/2<pp—Vss 

dVBS 


975 

QVbs 


■\/25P p — Vps  h  DFLTA 


Trg5/ 

ACoxW 


Cox 


(3.24) 


NFS  is  introduced  as  a  parameter  in  the  evaluation  of  Vq^.  Its  presence  is  neces¬ 
sary  to  invoke  the  weak  inversion  feature.  NFS  is  not  related  to  the  physical  nature  of 
the  subthreshold  conduction,  which  is  not  determined  by  the  fast  surface  state  density. 
It  is  a  curve-fitting  parameter  v;hich  can  be  extracted  from  measurements. 

In  the  above  equation  it  is  important  to  notice  that  the  quantity  y/Ztfip-Vps  which 
is  a  multiplicative  factor  of  the  depletion-cheirge  formulation  assumes  that  Vbs  is  nega¬ 
tive.  A  more  exact  way  to  express  Eq.  (3.24)  is  to  replace  the  square  root  by  a  F(  >55). 
In  the  new  M0S2  model,  this  function  is  made  to  decay  as3rmptotically  towards  zero 
when  Vbs  becomes  positive  rather  than  clamp  it  to  zero  as  v/-as  the  case  in  the 
SPICE2.D  and  .E  versions; 


The  weak  inversion  current  equation  for  Vcs<Von  is 


(3.25) 


Ids  =  P^^VoN-yBiN- 


V^DS 


Vr. 


PS 


-^/S'^S^^PF-ysS-^  ^£.5) 


3_  3. 

-{Z:p^-Vss)^ 


xe 


■^Vcs-Von) 


(3.26) 

This  equation  insures  the  continuity  of  the  current  at  Vqn  but  not  the  continuity  of 
the  first  derivative  in  contrast  to  the  previous  formulation  [2].  This  alternative  has 
been  chosen  since  the  earlier  expression  suffered  from  introducing  a  negative  resis¬ 
tance  in  the  Ids~Vds>Vgs  =  constant  characteristics. 


/ 
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In  order  to  obtain  the  usual  aspect  of  the  characteristics  in  the  los-Vjjs  plane  a 
saturation  voltage  is  defined  at  which  replaces  Vos  in  Eq.  (3.26)  when  the 

latter  becomes  larger  than  the  former.  The  role  of  the  exponential  factor  is  then  to 
scale  this  characteristic. 


3.5.  Temperature  Dependence 

The  M0S2  model  readjusts  all  temperature-dependent  variables  in  the  drain- 
source  current  equation.  These  include 

-  proportionality  with  T  of  the  Fermi  potential  <pF\ 

-  variation  with  temperature  of  the  energy  gap; 

-  temperature  dependence  of  the  built-in  voltage  PB  of  the  drain  emd  source  junctions; 

-  mobility  variation  as  T~^^\ 

-  temperature  variation  of  the  reverse  current  of  the  diffused  junctions. 


3.6.  Charge-Oriented  Model  for  HOSFETT  Cs^kacitances 

Meyer's  capaciteince  model  [13]  which  has  been  used  in  earlier  versions  of  SP1CE2 
suffers  from  a  fundamental  weakness:  it  does  not  conserve  chetrge.  An  improved  for¬ 
mulation  of  the  charge-oriented  model  proposed  by  Ward  and  Dutton  [14]  is  used  in 
SPICE2.G.  The  new  model  is  based  on  the  actual  distribution  of  charge  in  the  MOS 
structure  and  its  conservation: 

Qchan  -  Qd'^Qs  -  ~{Qg^Qb)  (3.27) 

The  current  fiowing  at  any  terminal  of  a  region  can  be  related  to  the  charge  contained 
in  that  region,  i.e.. 

. 

dt 


is  = 


dt 


(3.28) 


i-s+^D  - 


d(<?5  +  9x>) 
dt 


The  equivalent  terms  of  the  modified  nodal  admittance  matrix  are  found  by  using  the 
numerical  implicit  integration  algorithm  already  built  into  the  circuit  simulator  for  the 
transient  anadysis  [15], 
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(3.29) 


fidt  =  Q{t,)-Qito) 
^0 


The  charges  at  any  time  point  are  assumed  a  function  of  only  the  terminal  voltages  at 
that,  same  time  point.  Due  to  the  complexity  involved  in  finding  the  derivatives,  a  sim¬ 
ple  formulation  is  used  for  the  charges  [16]  which  is  adjusted  for  continuity  at  the  tran¬ 
sitions  between  regions  of  operation.  The  exact  formulations  are  the  subject  of  a 
future  write-up.  If  the  trapezoidal  integration  formula  is  used.  Eq.  (3,29)  can  be  rewrit¬ 
ten  as  follows 


IKVi+W  -  -9vo)+S 

X 


iVx\-yx°i) 


(3,30) 


where  x  stands  for  GB,  DB  and  SB,  and  y  for  G,  B,  D  and  S,  The  subscripts  indicate  the 
time  point  and  superscripts  the  iteration.  Thus  quantities  at  time  point  ”0"  are  known 
(previous  time  point),  quantities  at  time  point  ”1''  iteration  "0"  are  also  known  from  the 
solution  of  the  last  iteration  at  the  present  time  point.  Time.  "1”,  iteration  ”1"  is  the 
current  iteration  which,  by  substituting  the  above  equation  into  the  modified  nodal 
system,  has  V}i  as  a  solution. 

Defining  capacitances  for  the  small-s’gnal  analysis 


- 

•»*  "  dV^ 


(3.31) 


one  notes  that  between  each  pair  of  nodes  there  are  two  capacitances,  which  in  general 
are  different. 

Eig.  3.5,  plots  the  voltage  dependence  of  the  various  terminal  capacitances  associ¬ 
ated  with  the  Qc  and  Qg  charges.  Only  six  capacitances  oat  of  a  total  of  twelve  defined 
by  Eq.  (3.31)  are  independent  due  to  the  charge-conservation  principle  expressed  by 
Eq.  (3.27),  Qs  and  Qg  share  equal  parts  of  the  channel  charge  Qchan  in  the  linear 
region.  In  saturation  Qg  gets  only  XQC  x  Qchan  where  XQC  is  a  model  input  parameter 
which  defaults  to  zero;  Qs  results  from  Eq.  (3.27). 


1 
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4.  THi:  M0S3  MODEL;  A  SEMI^affiPIRICAL  MODEL  FOR  SMALL  GEOMETRY  MOSFET'S 


A  small-geometry  MOSFET,  defined  as  a  transistor  with.  L<Z}Mni,  and  W^Zjjjn  is 
characterized  by  the  following  features  : 

-  thres.hold-voltage  sensitivity  to  the  length  and  the  width  of  the  device  due  to  the  two- 
dimensional  nature  of  potential  distribution; 

-  threshold-voltage  sensitivity  to  the  drain  voltage  due  to  the  drain  induced  barrier 
lowering; 

-  relaxesd  transition  between  lineeir  and  satur^ation  regions,  and  lowered  saturation  vol¬ 
tage  and  saturation  current  due  to  the  velocity  saturation  o!  hot  electrons. 

The  M0S3  model  has  been  developed  to  address  the  above  features  euid  the  computa¬ 
tional  efficiency.  The  companion  capacitance  model,  which  conserves  charge  [14],  is 
also  derived  and  implemented. 

Instead  of  repeating  the  physical  explanations  already  presented  in  the  previous 
section,  the  M0S3  model  is  described  in  simple  terms  with  emphasis  on  the  different 
approaches  between  the  MOSS  and  MOS3  models,  followed  by  the  parameter  extraction 
methodology  of  those  newly  introduced  pexameters  and  the  performance  comparision 
between  M0S2  and  M0S3. 

All  the  parasitic  resistances,  overlap  and  junction  capacitances  are  described  by 
the  same  parameters  as  in  M0S2.  The  model  parameters  relating  to  the  intrinsic  MOS¬ 
FET  are  designed  to  be  compatible  with  those  of  M032.  The  rule  of  '’Electrical”  and 
"Processing"  parameters  also  holds  for  M0S3.  The  four  new  parameters  which  are 
specific  to  M0S3  are  described  in  Sec.  2.  ETA,  DELTA  and  KAPPA  are  dimensionless 
"Electrical”  pareimeters,  of  typical  values  close  to  1.  THETA  is  an  "Electrical"  parame¬ 
ter  of  unit  1/VoIt.  Its  typical  value  is  approximatedly  0. 1, 


4.1.  Threshold  Voltage 

The  relationship  between  Vtu,  the  threshold  voltage,  and  the  substrate  bias  devi¬ 
ates  from  the  first-order  square  root  dependence  as  the  device  dimensions  are  scaled 
dovm.  Tills  is  mainly  because  the  built-in  potential  depends  logarithmically  on  the  sub¬ 
strate  doping  which  is  scaled  up  proportionally  with  the  dimensions.  This  feature  is 
modeled  by  decoupling  the  effects  of  channeT length,  channel  width  and  static  feedback 
from  drain  to  gate.  Each  of  these  features  is  enabled  by  specifjdng  the  parameters  XJ, 
DELTA  and  ETA,  respectively. 

The  short-channel  effect,  the  decrease  of  as  the  channel  length  is  scaled 
down,  is  modelled  with  Dang's  [17]  modified  trapezoidal  approach,  reckoning  with  the 
influence  of  cylindrical  field  distribution  (see  Fig.  4.1). 
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Ilie  narrovA  channel  effect,  the  increase  of  Vj-f/  as  the  channel  width  is  scaled  down, 
is  modelled  by  taking  into  account  the  extra  bulk  charge  at  the  edge  of  the  channel  eis 
described  in.  the  previous  section.  To  include  other  edge  effects  such  as  the  existence 
of  field  implant  and  non-planarity  due  to  the  LOCOS  process,  the  empirical  parameter 
DELT^.  is  introduced. 

The  static  feedback  effect  can  be  explained  as  a  consequence  of  the  drain  induced 
barrier  lowering  [16],  A  term,  linearly  proportional  to  Vps,  is  included  in  the 
expresision.  The  constant  of  this  term  is  inversely  proportional  to  the  oxide  capaci¬ 
tance  and  the  cube  of  the  channel  length  [19].  The  parameter  ETA  is  introduced  to 
allow  more  flexibility. 

Ad.  expression  for  which  sums  up  the  above  features  is  formulated  as: 


Vjji  -  Vfg  H-  —  a  Vqs  +  yFs  +  -^v  {~9f~''3s) 


(4.1) 


where 


a  =  caefficwnt  of  static  feedback 


=  ETA 


n 

CqxL^ 


n  :=  empiricaJ.  constant 
■■=  8.15x10-2®  (/"xm) 

Fs  cQTTecticn  factor  of  short  channel  effect 


XJ 

LD^Wc 

i  Wp  /XJ 

*  LD 

L 

XJ  V  ^ 

[l  +  Wp/XJ 

XJ 

(4,2) 


(4.3) 


(4.4) 


LD  is  the  lateral  diffusion  length,  Wp  is  the  depletion  layer  width  of  a  plane  jvmctioa  Wc 
is  the  depletion  layer  width  of  a  cylindrical  junction  and 


— -,io  +  c£,— , 


(4.5) 


do,  dj  and  dg  are  empirical  constants  of  vjilues 


do  =  0.0631353 
di  =  0.6013292 
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dz  =  -0.01110777 

and 

Fn  =  correction  factor  of  narroiu  chmvnel  effect 


=  DELTA 


^si 

ZCqxW 


Xj)  -  coefficient  of  depletion  layer  vnjdth 


2g5J 

qNSUB 


(4.6) 


(4.7) 


4.2.  Basic  Drain  Current  Equation 

The  drain  current  can  be  expressed  as  : 


TDS=^^[vGS-VrH{y)\dVy 

^JOS 

=  P  f  [Vcs  ■'VrH-(l^FB)Vy]dVy 
0 

where 

^  =  ^r^EFF  ^OX 

Fb  -  Taylor  series  expansion  coefficient  of  bulk  charge 


7^s 

A'y/^p—VBS 


+  Fj/ 


(4.8) 


(4.9) 

(4.10) 


After  the  integration,  the  following  equation  is  obtained  [30]  : 

Vos  (4.11) 

With  such  simplicity,  an  explicit  saturation  voltage  expression  is  guaranteed.  Based  on 
the  assumption  that  the  Vos  dependent  term  in  Vth  represents  the  average  drain  vol¬ 
tage  influence  on  the  cheinnel  potential,  it  is  treated  as  a  constant  throughout  the' 
integration. 


hs  -  ^ 


Vos  ~  Vth - 
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4.n.  Surface  Mobility  Modulation  by  Gats  Voltage 

The  simplest  empirical  equation  is  used  to  enhance  the  computational  speed  : 


= 


_ UO _ 

1  +  THETA{V(s-Vth) 


(4.12) 


4.4  .  Velocity  Saturation  of  Hot  Electrons 

The  saturation  of  hot  electron  velocity  lowers  the  conduction  current  in  the  linear 
region.  Thus  relaxes  the  transition  between  the  linear  and  the  saturation  regions.  In 
the  linear  region  this  effect  is  modelled  by  the  conventional  hyperbolic  equation  [21]  : 


P-EFF  - 


1  + 


niAXxL 


'DS 


(4.13) 


If  tlae  parameter  VMAX  is  not  specified  by  the  user,  is  set  to  fis  and  this  effect  is 
not  modelled. 


4.5.  Saturation  Voltage 

The  saturation  voltage  of  a  short-channel  device  is  the  drain  voltage  at  which  the 
carriers  asymptotically  reach  the  maximum  velocity  at  the  dredn  limit  (see  Fig.  4.2)  : 

IBs  =  Qui^rain)  VMAX  (4.14) 


After  the  substitution  of  the  expressions  of  Qr^  and  iBs.  the  above  equation  becomes  : 


2 


'^DSAT 


~  ^TH  —(l+^s)  ^DSAT 


VMAX 


(4.15) 


From  the  above  equation  Vos  at  is  solved, 


'^DSAT  - 


l^GS  ~  %/ 
1  +  Ffi 


yMAXy-L 

f^s 


Vi- 


1 

'  VMAX^L^ 

l  +  F^ 

i^s 

(4.16) 


^DSAT  'Will  depend  on  Vqs  if  the  parameter  ETA  is  not  zero.  F  parameter  VMAX  is  absent 
from,  the  input,  the  saturation  voltage  will  be  determined  by  the  maximum  of  the  drain 
current  equation,  corresponding  to  the  situation  of  channel  pinch-off  : 
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(4.17) 


^DSAT  - 


1^C5  ~  yfH 
1  +Fg 


4.6.  Channel  Length  Modulation 

As  Vjxs  becomes  greater  than  YdsaTi  tbe  point  where  the  carrier  velocity  begins  to 
saturate  moves  towards  the  source.  The  channel  length  reduction,  AL,  is  formulated 
based  on  Baum’s  theory  [7]  : 


At 


+  KAPPAxXS  {Vds-Vdsat)  - 


XSEp 

2 


where 


Ep  -  lateral,  field  at  channel  pinch —off  point 


(4.18) 


(4.19) 


hsA7 


hsAT  and  Gbsat  are  the  drain  current  and  the  drain  conductance  at  saturation  respec¬ 
tively. 

It  should  be  noted  that  the  point  where  the  velocity  begins  to  saturate  differs  from 
the  point  v/here  the  inversion  charge  goes  to  zero,  i.e,.  the  channel  pinch-off  point. 
ycsAT  is  the  voltage  at  velocity  saturation,  while  Ep  is  the  lateral  field  at  channel 
pinch-off.  Thus  the  voltage  across  the  depleted  surface,  of  length  AZ,  should  be  less 
than  {Vbs  -  Vdsat)-  The  parameter  I^AI^PA  is  introduced  to  account  for  this  effect.  As 
AZ,  appraoches  I.  the  same  scheme  as  described  in  Sec.  3.3.1  is  used  to  prevent  a 
negative  effective  channel  length. 


4.7.  Weak  Inversion  Conduction 

The  equations  in  the  weak  inversion  region  are  the  same  as  those  used  in  M0S2, 
which  provide  both  the  continuity  of  the  drain  current,  the  proper  bias  dependence 
eind  sufficient  computaticnal  efficiency.  The  same  parameter  NFS  is  used  to  turn  on 
the  vreak  inversion  feature.  For  details  see  Sec.  3.4. 
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4.8.  Temperature  Dependence 

All  the  quantities  referred  to  in  the  model  equations  have  the  same  temperature 
dependence  as  in  M0S2.  The  four  newly  introduced  parameters  are  all  empirical.  No 
temperature  dependence  has  been  developed  for  them  as  yet. 


4.9.  Caqpacitance  Model  vrith  Charge  Conservation 

The  total  amount  of  charge  residirig  on  the  gate,  Qo,  can  be  evaluated  as 


Qa  =  W  fqg{y)dy 


(4.20) 


where 

qg  =  gale  charge  per  unit  area 
=  Cqx  Ifgs  -  {Vf^  +  2cpF-a  Vos)  - 
q^  =  channel  charge  per  unit  area 


(4.21) 


(4.22) 


-  ~Oox 


VgS  -  Vth  -  {'^■hFg)  Vy 


After  carrying  out  the  integration, 


Qa  =  WxLxCax 
where 

Fi  =  Vos  -  V'rn 


Vcs  ~  ( Vjro  -t-  2<pp  —  a  Vos)  ~ 


T^5  1  +  Fb 


12  Fi 


Vh 


DS 


^^Fb  . 


DS 


Similarly,  one  can  get  the  total  bulk  charge  Qb  as 
Qb  ~ 


Fb 


[y jTs  +  Fb  (2iPb-Vbs)  +  -^Vos 


Fb(1  +  Fb) 
12  Fj 


V$s 


(4.23) 


(4.24) 


(4.25) 


The  total  channel  charge  Qcbmi  results  from  the  charge  conservation  principle, 


Qohan  -  -{Qg-^Qb)  (4.26) 

It  is  equally  divided  into  the  source  and  the  drain  charges  in  the  linear  region  of 
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operation.  In  saturation  Qo  gets  a  share  XQC  of  Qchan- 

ilach  of  these  charges  has  three  associated  derivatives  ’which  are  the  capacitive 
elements  in  the  circuit  model.  It  should  be  noted  that  only  six  of  these  twelve  capaci¬ 
tive  components  are  independent. 


4. 10.  Parameter  fbriractioa  for  HOS3 

Tlie  characteristics  of  small  geometry  MOSFET's  are  sensitive  to  the  device  dimen¬ 
sions.  The  systematic  parameter  extraction  requires  measurements  of  several  devices 
with  different  combinations  of  length  and  width,  eind  the  manipulation  of  the  collected 
data.  The  fundamental  parameters  VTO.  NSUB,  UO  and/or  GMIMA  and  KP.  which  are 
common  to  both  M0S2  and  M0S3,  should  be  extracted  from  a  lOiig  and  wide  device,  as 
described  in  next  section.  Only  the  extraction  methodology  of  parameters  relating  to 
the  second-order  effects  is  described  in  the  following  : 


4. 10. 1.  THErA.  Mobility  Modulation  Parameter 

The  effective  surface  mobility  is  the  slope  of  the  [(Ix/ij5)/(  versus  Fns 

plot  in  I'.he  linear  region,  at  constant  Vj)s  and  Vgs,  as  indicated  by  the  current  equation, 
(4.11).  Rearrange  Eq.  (4.12)  as  ; 

-^■=  1  +  THETAiVas-VrH)  (4.27) 

Clearly,  parameter  THETA  is  the  slope  of  the  {UO/^s)  versus  Vgs  plot.  In  order  to  iso¬ 
late  the  gate  voltage  influence  from  the  hot  electron  effect,  the  measurements  of  Igs 
versus  Fas'  must  be  carried  out  at  low  Vgs,  say  0.2  volts. 


4.10.2.  ISTA,  Static  Feedback  Parameter 

The  coefficient  a,  the  slope  of  the  Vjg  versus  plot  at  constant  Vgs  as  defined  in 
Eq.  (4.2)  represents  the  magnitude  of  the  static  feedback  effect.  The  determination  of 
a  involves  the  measurement  of  the  threshold  voltage  of  small  geometry  MOSFET's  at 
interrneciate  Fi,s  level.  The  should  be  determined  from  the  Igs  versus  Vgs  curves 
at  different  at  low  current  levels  (see  Fig.  4.3). 


vol^eig*  o‘F'P««‘b 


d«ll;.ei  VTH 


VGS-VFB  <volt.> 
- yos-i.  0 _ j/Ds-a.  0 


Other  threshold  voltage  determination  algorithms,  based  on  either  the  simple 
saturation  current  equation,  which  is  proportional  to  the  square  of  {Vqs-Vth)  or  the 
approximately  linear  current  equation  v^hich  is  proportional  to  705/2)1^5. 

does  not  hold  because  of  the  hot  electron  effect.  Parameter  ETA  is  related  to  a  by 


ETA=a 


Q 


(4.28) 


The  value  of  the  empirical  constant  fl.  see  Eq.  (4.3),  has  been  adjusted  such  that,  for 
most  practical  devices,  ETA  assumes  a  value  close  to  one. 


4. 10.3.  DELTA,  Width  Effect  on  Threshold  Voltage 

Tlie  coefficient  of  the  width  effect  on  Fj®  at  constant  relates  inversely  to  the 
channel  width  as  indicated  by  Eq.  (4.6).  DELTA  relates  to  F}f{Ztpjp~VjBs)ty ,  the  slope  of 
Vm  versus  1  /W  plots,  as  : 


DELTA  -  slope 


^  Oar _ 1 

rrssi  1^5) 


Its  value  varies  from  process  to  process. 


(4.29) 


4.10.4.  VMAX.  limited  Velocity  of  Hot  ESectrons 

The'  effective  mobility  at  intermediate  Vj)s  level  deviates  from  the  value  at  low 
Ps-  Thc!  parameter  VMAX  can  be  estimated  as  follows  :  measvure  [j)s  versus  Vac  at 
different,  fixed  extract  the  effective  mobility,  psFT’  from  the  collected  data,  and 

finally  plot  (l/MaPF— 1/Ms)  versus  Viu/L.  The  slope  of  this  plot  is  the  inverse  of  VMAX 
as  indicated  by  the  rearrangement  of  Eq.  (4.13)  : 


1  I  _  Fps 
pEFf  Ps  VMAX  X  L 


(4.30) 


4. 10.5.  K4Pi^A.  Correlation.  Coefficient  of  the  Drain  field  in  Saturation 

The  extraction  of  parameter  KAPPA  from  measurement  is  not  as  straight  forward 
as  for  the  other  parameters.  First,  one  has  to  estimate  the  approximate  saturation 
voltage  ar.d  saturation  current  using  the  VMAX  obtained,  as  above.  Then  the  channel 
length  reduction  can  be  estimated  as  : 

M  =  (4.31) 

hSAT 

The  slope  cf  the  versus  plot  is  approximately  KAPPAxXS, 
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4. 11.  Performance  Comparison  betireen  MQS2  and  M0S3 

Though  M0S3  and  MOSS  share  most  o'  the  model  parameters,  different  values  for 
the  same  parameter  must  be  used  to  produce  approximately  the  same  characteristics. 
For  example,  the  parameter  VMAX  has  no  effect  on  the  linear  region  characteristics  of 
MOSS,  while  it  lowers  the  effective  mobility  of  MOSS.  In  order  to  get  a  match,  a  lower 
value  for  UO  in  MOSS  is  required. 

The  Vbs  dependence  of  the  basic  drain  current  equation  in  MOSS  is  an  approxima¬ 
tion  of  that  in  MOSS.  With  all  other  parameters  being  properly  adjusted  for  a  specific 
device,  say  L=5/i,m,JV5C/F=4xl0*®cm~^andrQY=0.1/im.  this  approximation  only  intro¬ 
duces  about  3%  deviation  in  the  case  of  Vqs  well  above  Vf^-  The  deviation  increases  to 
S0%  when  the  gate  voltage  is  only  0.5  volts  above  Vth-  However,  In  the  latter  case,  the 
device  is  already  approaching  subthreshold  conduction  and  the  basic  strong  inversion 
theory  begins  to  fail. 

The  benchmark  runs  with  SPICES. G  indicate  that,  with  compatible  input  parame¬ 
ters,  the  MOS3  model  is  up  to  40%  faster  in  the  model  computation  time  as  compared 
with  MOSS,  depending  upon  the  nature  of  the  circuit  and  its  operation. 
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5.  MODEL  PARAMETER  EXTRACJION 

Naturally  the  accuracy  of  the  results  of  a  MOS  IC  simulation  depends  heavily  on 
the  values  of  the  model  input  parameters:  The  I-V  cheu’acteristics  predicted  by  the 
equations  of  the  model  must  match  very  closely  the  measured  characteristics.  The 
coraplexity  of  the  equations  in  Sec.  3  is  such  that  some  computer  aids  are  required  to 
obtain  the  program  value  of  the  device  p^r^eters.  The  best  results  are  obtained  if  a 
characterization  system  is  set  up  [22]  whi,ch  has  as  core  a  desk-top  or  a  mini  computer 
which  controls  the  measuring  instrument  and  does  the  necessary  computations  for 
fincLing  the  model  parameters.  The  equaUons  implemented  in  the  parameter  extrac¬ 
tion  program  must  be  identical  to  thps|e  used  in  the  circuit  simulator.  A  data  acquisi¬ 
tion  software  designed  for  a  certeun  model  ppoviles  accurate  input  parameters  only  for 
that  particular  model.  When  the  model  is  changed,  the  parameter  extraction  software 
has  to  be  changed  as  well. 

Another  tool  necessary  for  parameter  extraction  is  a  test  chip  which  should  have 
the  following  features: 

-  a  short-  and  a  long  channel  MOSFET; 

-  a  Jiarrow-  and  a  wide  channel  device; 

-  dijJused  resistances; 

-  pcily-silicon  resistances; 

-  a  I’ectangular  and  a  meander  form  junction  capacitor; 

-  a  thin  and  a  field  oxide  capacitor. 

A  desired  feature  of  the  parameter  extraction  software  is  to  find  the  spread  of  the 
parameters  for  a  number  of  wafer  runs  and  to  group  the  parameters  according  to  a 
certcun  criterion,  e.g.,  speed  or  power  consumption,  into  a  worst,  best  and  nominal 
caS'S  input  file  [23]. 

It  is  the  purpose  of  this  section  to  illustrate  the  methodology  and  the  minimum 
hardware  software  to  extract  the  main  input  parameters  for  the  M0S2  model.  The 
reported  work  hcis  been  performed  using  a  Tektronix  4051  desktop  calculator,  Tek¬ 
tronix  576  curve  tracer,  a  C-V  measuring  instrument  and  a  multimeter. 


5. 1.  VTO  and  NSUB  Measurement 

The  measurement  is  taken  on  a  saturated  transistor  connected  to  a  curve  tracer. 
Sim:e  at  this  starting  point  few  parameters  are  known,  the  saturated  characteristic  is 
approximated  by  a  parabola  rather  them  using  the  exact  current  equation  given  in 
Sec  3. 
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where  Qj  is  the  ordinate  at  the  origin  of  the  C-V  plot  and  Ai  is  the  area  (known)  of  the 
rectangular  capacitor.  The  measured  resvilts  indicate  that  MJ  is  very  close  to  0.5. 

CJS¥  and  MJSW  can  be  found  using  a  linear  regression  through  the  {Vr.Cj^  data 
points  of  the  meander-form  diffusioa  Eq.  (5.4)  can  be  rewritten  eis 


log  Ck  =  log  (PxCJSW)  -  MJSWxlog 


(5.5) 


where 


CJ 


1+ 


PB 


ITT 


(5.6) 


The  data  processing  is  performed  by  PROGRAMl  (see  AppendiK  1)  which  in 
response  to  the  point  pairs  (Vr.Cjs)  and  other  constants  displays  the  result  of  the 
interpolation,  i.e.,  CJSW  and  MJSW. 

If  MJ  =  0.5  is  not  adequate,  PROGRAMl  can  be  easily  modified  to  perform  the  same 
linear  regression  as  Eq.  (5.5)  for  the  large  rectangular  capacitor  only  to  find  MJ. 


5.3.  Mobility  Parameter  Extracti  on 


Finding  the  parameters  which  describe  the  mobility  variation  with  the  surface 
field  is  very  important  for  am  accurate  simulation  of  the  I-V  characteristic.  UO,  UCRIT 
and  UEXP  which  appear  in  the  mobility  degradation  formula  of  M0S2  (see  Eq.  (2.7)), 
must  be  found. 

The  experimental  setup  allovrs  independent  variedion  of  Vqs  and  Fes'  on.  the  curve 
tracer  and  to  extract  data  points  from  the  linear-  (very  low  Vjc)  and  triode  region  up 
to  VdsaT' 

A  straight  line  can  be  interpolated  through  these  data  points  by  rearranging 
Eq.  (2.7) 


log 


=  UEXP. 


log  (-^C/CR/T)-log  {Vcs-Vth-UTF.A  x  V^s) 
'-'OX 


(5.7) 


where 
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(5.1) 


The  threshold  Vm  is  taken  as  the  intercept  on  the  abscissa  of  the  above  lineeir  regres¬ 
sion  through  the  coordinate  pair  (V^,  Vos).  Several  Vth  are  calculated;  one  for  each 

^5. 


VfH  =  WO + GAMMA{  ■sJVsb-^PHI-'JTW) 


(5.2) 


where  PHI  =  Zfp.  From  another  straight  line  interpolated  in  the  plane 
( V  +PHI,  Vj’ff') ,  the  slope  gives  GAMMA  and 


NSUB- 


GAMMA^ 

2965/ 


(5,3) 


All  the  above  extrapolations  are  performed  by  PROGItAMl  written  in  BASIC  for  the 
Tek  4051  the  listing  and  operating  instructions  of  which  are  contained  in  Appendix  1. 
The  user  has  mainly  to  input  only  the  {Ids  ,  Pa?)  point  pairs  and  the  1^5  stepping  and 
the  pr  ogram  will  display  VTO,  GAMMA  and  NSUB. 

These  measurements  should  be  performed  on  a  long  and  wide  channel  MOSFET 
since  the  above  values  will  then  be  corrected  by  the  program  for  small-size  effects,  For 
more  accuracy  the  extraction  program  could  perform  one  more  iteration  once  the  first 
estimate  has  been  found  for  VTO  and  NSUB  and  use  the  exact  current  equation  as 
implemented  in  M0S2. 


5.2.  Capacitance  Measurements 

As  has  been  already  mentioned  in  Sec,  2  the  jimction  capacitance  is  separated  into 
a  bottom  and  a  sidewall  capacitance 


CJ 


1+ 


Yl. 

PB 


W 


+P 


1  + 


CfSW 


PB 


(5.4) 


where  A  represents  the  plane  area  (on  layout),  P  the  perimeter  and  Vjf  the  reverse  bias 
applied  to  the  junction. 

CJ  can  be  found  first  by  assuming  that  the  sidewall  capacitance  of  the  large  rec¬ 
tangular  capacitor  is  negligible.  From  a  C-V  plot,  one  gets 
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(5.B) 


/^5  _  ^ i 

UO  ~  KPo 


with 

KPi  -  conduction  factor  for  different  V^y 

KPo  -  maximum  value  of  conduction  factor 

UEXP  and  UCRIT  are  obtained  from  the  above  linear  regression  as  the  slope  and  y 
intercept  respectively  by  plotting  log{KPi  /  KPq)  on  the  y-  and  lQgiVcs—Vjj/  —  UTPAxVi}s) 
on  the  X  axis. 

The  data  processing  is  performed  again  by  PROGRAMl  by  specifying  Ids) 

points  where  for  different  Vqs  curves  different  maximum  Vus  values  are  entered  in 
order  to  cover  each  characteristic  up  to  V^sat-  PROGRAMl  computes  UO.  UCRIT  and 
UEXP  for  each  Vj)s  range  and  at  the  end  displays  an  average  for  each  of  the  three 
values  over  the  whole  interval. 


5.4.  i^proximation  of  VMAX 

A  second  program  written  in  BASIC  tor  the  Tek  4051,  PROGRAMS  (see  Appendix  2) 
helps  with  its  two  options  to  evaluate  VMAX.  The  first  alternative  is  to  plot  V^sat  as  a 
function  of  the  effective  channel  length  L  with  Ucs  as  a  parameter.  For  each  VMAX,  one 
gets  a  different  set  of  curves.  Fig.  5.1.  shows  such  a  sample  plot. 

The  second  option  plots  Idsat  as  a  function  of  VMAX  with  as  a  parameter.  A 
different  set  of  curves  results  for  each  L.  A  sample  plot  for  this  option  is  shown  in 
Fig.  5.2. 

Starting  from  the  measured  characteristics  one  defines  Vdsat  as  the  voltage  start¬ 
ing  from  which  the  I-V  curve  can  be  approximated  by  a  straight  line  and  the  current 
corresponding  to  this  point  is  Idsat-  i^ith  these  values  one  can  get  a  good  guess  for 
VMAX  from  plots  similar  to  those  in  Fig.  5.1.  and  5.2.  It  is  emphasized  that  PROGIiAM2 
uses  the  same  equations  as  those  implemented  in  SPiCEE.G  for  Vdsat  (VfvIAX). 
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6.  SIMULATION  EXAMPLE 


It  is  the  purpose  of  this  example  to  show  the  ease  with  which  a  section  of  a  CMOS 
LS  I  circuit  is  accurately  simulated  starting  from  the  layout  and  using  the  new  device 
cerd  format  and  model  input  parameter.  The  circuit  of  Fig.  6.1  represents  a  short  path 
of  a  MOS/LSI  circuit.  The  output  signal  CS  generated  from  CLK  can  be  measured  at  an 
01  tput  pad  and  compared  with  the  simulated  waveform.  The  input  deck  for  SPICE2.G  is 
gi  ^en  in  Appendix  3. 

The  measured  circuit  chips  were/located  on  the  same  wafer  with  the  test  patterns 
fri>m  which  the  parameter  extractionihas  been  performed.  For  an  accurate  simulation 
ai  the  parasitics  of  the  layout  must  be  Included;  /?i  represents  the  series  resistance  of 
the  input  pad  protection  while  i?2  and  R3  represent  crossunder  diffusions.  Each  C 
indudes  all  the  capacitive  effects  associated  with  that  node,  e.g.,  Ci  is  contributed  by 
Al-  to-neld  oxide,  Poly  to-field,  N+  and  P+  junctions  and  thin  oxide,  where  for  each  the 
accurate  overlap  area  must  be  evaluated.  Each  of  the  gates  has  been  modeled  as  a 
subcircuit.  A  load  capacitance  of  lOOpF  has  been  used;  the  rest  is  attributed  to  the 
probe  and  parasitics.  Fig.  8.2.  shows  a  sample  simulated  output  of  the  waveforms. 

An  interesting  result  has  been  obtained  by  performing  the  simulation  with  and 
without  sidewall  capacitance  and  with  and  without  scattering-limited  drift-velocity 
efiects  (YMAX).  The  input  file  given  in  Appendix  3  is  complete  in  the  sense  that  it 
in  dudes  series  resistance,  sidewall  capacitance  and  drift  velocity  satioration.  The  dev- 
ictys  which  have  no  perimeter  specified  had  a  negligible  sidewall  compared  to  the  area. 

The  results  of  the  different  simulations  and  the  measured  values  are  summarized 
in  Table  I.  It  can  be  noticed  that  the  above  two  effects  contribute  a  15  to  20%  accuracy 
irrprovement  and  the  simulated  results  of  this  circuit  containing  22  MOSFET's  is  within 
5  to  10%  of  the  results  for  the  actual  circuit. 
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FIG.  6.2. 


TABLE  I 


SIMULATED 

PARAMETER 

MEASURED 

W/0  Sidewall 

W/  Sid 

ewa  1 1 

W/C  VMAX 

w;  VMAX 

W/0  VMAX 

W/  VMAX 

Delay 

■1 

64  -  75 

49 

51 

60 

62 

T„  (Nj) 

90  -  105 

70 

79 

81 

91 

Fall  Time 

C|_  -  20pF 

32  -  36 

22 

24 

24 

26 

Tp  (Nj) 

=  120pF 

P5  -  120 

90 

100 

94 

102 

i^pendix  1 

User  Instructions  and  Listing  of  PR0GRA.M1 

Load  the  program  from  the  tape  cartridge  (Tektronix  4051)  by  typing: 

FILE  n  CR 
OLD  CR 

^nX3,  GAMMA,  and  NSUB  Intei^lation 
Calculator  prompts; 

TYPE  OI”  UN  REGRESSION,  DEVICE  TYPE. 

NO.  OF  POINT  PAIRS,  NO.  OF  CURIAS 
Ansv,er 


GAZMMA 

CR 

N  (or  P) 

CR 

no.  of  ilDS’^as)  pairs  for  each  V^h. 

no.  of  different  VSB  values 

CR 

Prompt: 

PHI,  TOX,  VSB  STEP 

Enter  appropriate  signed  values;  although  NSUB  is  not  yet  known  specify  guess  for  PHI 
Prompt: 

lu,  VSB  =  value,  VGS  =  ,  IDS(UA)  = 

where  n  specifies  the  coordinate  pair  number,  value  the  1^5  value  and  the  user  enters: 

VGS  valvie,  IDS  value  in  /jA  CR 

This  is  repeated  as  many  times  as  necessary  for  entering  aU  data  as  specified  at  the 
begirning;  VGS  must  have  the  appropriate  sign. 

At  conclusion  the  program  will  print: 

WO  =  value,  GAMMA  =  value,  NSUB  =  value 


-45- 


Ca{>etcitance  Interpolation 


Answer  in  response  to  the  starting  prompt  (TYPE  OF  LIN  REGRESSION,  etc); 


CAPACITANCE 

CR 

N  (or  P) 

CR 

no.  of  (I^  ,  Cjs)  point  pairs. 

1 

CR 

Prompt: 

PB,  CJ  (PF/U-2),  AREA  (U-2).  PERIMETER 

where  PB  and  CJ  have  the  same  meaning  as  In  SP1CE2  and  AREA  and  PERIMETER  refer 
to  the  meander  form  capacitor. 

Enter  the  corresponding  values  separated  by  commas  and  CR  after  the  last  one. 

At  the  conclusion  the  calculator  vrill  display; 

CJSW  =  value.  MJSW  =  value 

The  program  assumes  CJ  knovm  and  interpolates  for  MJ=0,5;  it  can  be  easily  modified 
for  a  two  step  interpolation,  first  for  CJ  and  MJ  with  the  area  of  the  rectangular  capaci¬ 
tor,  and  then  for  CJSW  and  MJSW. 

HoMlity  Interpolation 

Answer  in  response  to  the  starting  prompt  (TYPE  OF  LIN  REGRESSION,  etc); 

MOBILITY  CR 

N  (or  P)  CR 

no.  of  (VgsJos)  point  pairs, 
no.  of  different  VDS  ranges  CR 

rhe  range  of  Vqs  has  the  meaning  that,  different  Vj)s  values  can  be  specified  for  each 
/Gj  =  constant  characteristic  in  order  to  cover  each  V^s  curve  up  to  the  saturation 
•oint  VjjsAT-  For  the  sample  data  shovm  in  Table  A.1  enter; 


MOBILITY 

CR 

P 

CR 

4,4 

CR 

ompt: 

VTO,  NSUB,  TOX,  UTRA.  W,  LEFF,  VSB 

respective  values  are  entered  separated  by  commas  and  a  CR  after  the  last  one. 
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LEFF  is  the  layout  channel  length  less  the  under-diffusions. 

Prompt: 

n.  VDS  =  ,  VGS  =  .  IDS(UA)  = 

where  n  is  the  current  number  of  data  point  set;  at  each  such  prompt  one  line  from 
Table  A.I  has  to  be  input.  After  all  data  for  one  Vds  range  have  been  input  the 
corresponding  UO.  UCRIT,  and  UEXP  are  displayed: 

YDS  =  value,  UO  =  value,  UCRIT  =  value,  UEXP  =  value 

This  operation  is  repeated  until  all  the  different  ranges  are  entered.  At  the  end  the 
program  will  print  an  average  value  for  the  three  mobility  factors  over  all  the  data 
points.  An  example  is  shown  in  Fig.  Al. 


TABLE  AI. 


B 

VGS 

IDS  (uA) 

11^ 

I^B 

I^B 

H 

-4.0 

30 

-0.4 

-2.5 

32 

-3.0 

40 

-3.5 

49 

-4.0 

57 

-1.0 

-3.0 

77 

-3.5 

102 

-4.0 

120 

-4.5 

145 

125 

180 

220 

-5.0 

255 
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90  INIT 

100  LET  G^=" GAMMA" 

110  LET  U$=''M0BILITV" 

120  LET  C$="CAPACITANCE“ 

130  PRINT  "TYPE  OF  LIN  REGRESSION,  OEUICE  TYPE," 

140  PRINT  "NO.  OF  POINT  PAIRS, NO.  OF  CURUES" 

150  INPUT  Tf,Q$,M,Ml 
IbO  T=1 

170  IF  Q$="N"  THEN  190 
ISO  T=-l 

190  DIM  K 1 < M ) , Ub < M > , X  < M > , Y  < M > , U4 ( M > 

200  D I M  U0< M 1 > , U 1 ( M 1 > , U2 C M 1 > , U9 < M 1 > , B2 CM 1 ) 

210  N9=1.45E+10 
220  N1=0 
230  N=0 
240  N1=HH1 
250  N=N+1 

200  IF  Tf=C$  THEN  1350 

270  IF  Tf=G$  THEN  920 

230  IF  T$=U$  THEN  310 

250  PRINT  "HOW  ABOUT  ANOTHER  READING  OF  THE  USERS 'MMANUAL,  CUTIE?" 
300  STOP  r 

310  IF  N>1  THEN  440 
320  IF  N1>1  THEN  440 

330  PRINT  "UTO,NSUB,TOX<U>,UTRA,M<U),LEFF<U>,USB" 

340  INPUT  U1,N0,T0,U3,M,L,U5 

350  T0=T0T1.0E“4 

360  W=WT1.0E-4 

370  L=L4c1.0E~4 

380  E2=1.04E-12 

390  EB=8.845E-14^:3.9 

400  F2=TT0.052TLOG(N0/N9) 

410  G0=SQR<3.2E“l9*E2tH0)/E0*T0 
420  U0=Ul'-F2-TTG8JrSQR<T*F2> 


PROGRAM  1  (continued) 


430  A0*T*<F2+U5) 

440  PRINT  N,"UDS='',“UGS=% ''IDS<UA)  =  “ 

450  INPUT  U4(N>,U6(N>,Y0 
460  A1=A0+TTU4<N) 

4?0  A2=SQR(A0>t3 
4S0  A3=SQR<Al)t3 
490  E!0=2/31G0t(A3-A2> 

500  1 0=  < U04F2+U4 ( N > /2> *04  < N ) 

510  U2=Ul+U3tU4(N> 

520  Il=Y0t.l.0E"D 

530  D0=U6 ( H >  t U4 ( H  > - 1 0-B0 

540  K1(H>=I1/D0 

550  IF  N=1  THEN  570 

560  IF  N>1  AND  K1(N><=K0  THEN  580 

570  K0=K1(H> 

5S0  IF  N<M  THEN  250 
590  S1=0 
600  S2=0 
610  S3=0 
620  S4=0 

630  FOR  J=l  TO  H 

640  K1<J>=LGT(K1(J>/K0) 

650  U6(J>=LGT(T.t(06<J>“Ul~U3.*U4(H>>) 

660  S1=S1+U6(J) 

670  S2=S2+UG(J)t2 
680  S3=S3+K1(J> 

690  S4=S4+U6<J>tKi(J) 

700  NEXT  J 

710  U0(N1>=K0*T0/E0*L/W 

720  N0= (83*82-5 ltS4)/<H*S2-Slt2) 

730  S0=(N*S4-S1*S3>/(N*S2-S112) 

740  U2<H1>=-S0 

750  UKNl  )  =  10t(-N0/S0>*3.9/<T0*11.7) 

760  PRINT  ••UDS=";U4(N>r'  U0=”  i  U0(NI )  i  -  UCRIT*"  J U1  <N1  >  J  UEXP=- J U2<N1  > 
770  IF  NHHl  THEN  230 


I 

cn 

ro 

I 


780 

790 

800 

810 

820 

830 

8^0 

850 

880 

870 

830 

890 

900 

910 

920 

930 

940 

950 

960 

980 

990 

1000 

1010 


U6=0 

U7=0 

US=0 

FOR  K*1  TO  Ml 
Ub=U6+U0<K> 
U7=U7HI1(K) 
U8=U8^-U2CK) 
NEXT  K 
U6=Ub/Ml 
U7=U7.''M1 
U8=U8/M1 
PRINT 


15"^^ 


•tt.t  AVERAGE  OUER  S 
U0=“;U6,  "UCRIT 

u  r  0 

THEN  1060 

AND  H1>1  THEN  1000 
"PHIjTOXiU.MiSB  STEP" 


PRINT 
GO  TO 
IF  N>1 
IF  N=1 
PRINT 

INPUT  F2,T0,D1 
T0=T0.1:1.0E-4 
U5=0 

GO  TO  1010 
U5=l)5+Dl 
S1=0 


PECIFIED  IHTERUAL  tU 
=";U7,"UEXP=";U8 


It 


1020  S2=0 


1030  S3=0 
1040  S4=0 

1050  B2<N1 )sSQR<T*(U5+F2) ) 
1060  PRINT  H» "USB-";U5i "UGS*" 
1070  INPUT  X0,Y0 
1030  Y2=SG!R<Y0T1.0E-6> 

1090  Sl=Sl+T:|fX0 
1100  S2=S2+X0t2 
1110  S3=£3+Y2 


1120  S4s:S4+TTX0*Y2 
1130  IF  N<M  THEN  250 


••IDS<UA>  =  " 


PROGRAM  1  (continued) 


1140  D2=N1:S2-Slt2 
1150  N0=<S3*S2“S1*S4)/D2 
1160  S0=<H1:S4-S11:S3>/D2 
1170  U9<N1)  =  -T1:NO.''S0 
1180  IF  Nl<ril  THEN  230 
1190  S1=0 
1200  S2=0 
1210  S3=0 
1220  S4=0 

1230  FOR  J=1  TO  Ml 
1240  Sl=Sl+B2(J> 

1250  S2=S2+B2(J>t2 
1260  S3=S3+TtU9(J> 

1270  S4=S4+TTU9<J>*B2(J) 

1280  NEXT  J 

1290  D2=»'m^S2-Slt2 

1300  H0=(S3t.S2-Sl*S4)/D2 

1310  G2=<Nir.S4-Slt.S3>/D2 

1320  ll5=(3.3tG2/T0>T2t8.854E“14/3.744E-18 

1330  PRINT  "UT0=";li9(l>,  ''GAMMh='';G2,'‘NSUB="jN5 

1340  GO  TO  1570 

1350  IF  H>1  THEN  1420 

1360  PRINT  ■'PB,Cfi0<PF/Ut2>,AREA<U12>,PERIMETER<U>'' 

1370  INPUT  F3,C0,AO,P0 

1380  S1=0 

1390  S2=0 

1400  S3=0 

1410  S4=0 

1420  PRINT  N, "URjCJCPF)" 

1430  INPUT  U1,C1 

1440  K=C1-A04CO/SOR<1+U1/F3) 

1450  X0=LGT<l-»-Ul/F3) 

1460  Y0=LGT(K) 

1470  S1=S1+X0 
1480  S2=S2+X0t2 


I  ISUUIMMI  i 


vuoiiL  iiiueu; 


1490  S3*S3+Y0 

1500  S4=S4+N01:Y0 

1510  IF  N<H  THEN  240 

1520  Net=(S3l;S2-SUS4)^'(Hl:S2“Slt2) 

1530  S0=<NtS4-S11:S3)/<NtS2-Slt2) 

1540  C2=10t(N0-6>/P0 

1550  I'1O=-S0 

1560  PRINT  “CP0i;F/M)=  “;C2,’‘M=  “IHO 
1570  END 


Appeadbc  2 

User  Instructions  and  Listing  of 

Use  the  same  procedure  as  in  Appendix  1  to  load  the  program  from  the  tape  cartridge. 
Prompt; 

SELECT  TYPE  OF  PLOT:  ENTER  1  FOR  VDSAT  =  F(L) 

2  FOR  IDSAT  =  F(VMAX) 

Enter: 

1  (or  2)  CR 

A  first  set  of  parameters  is  required  next; 

TYPE,  VTO,  NSUB,  TOX(U).  W(U),  LEFF(U),  VSB 
Enter  the  required  data  as  follows: 

N  (or  P)  CR 
values  CR 

The  values  are  input  separated  by  commas  eind  a  CR  at  the  end  of  the  line. 

Additional  parameters  are  required  next: 

UO.  UCRIT,  UEXP,  UTRA,  VMAX  (M/S) 

These  are  input  in  the  same  dimensions  as  used  inSPICE2.G. 

The  Umits  and  step  size  of  the  variable  (L  or  VMAX  depending  on  wliich  t3q>e  of  plot  has 
been  selected,  1  or  2  respectively)  are  requested  next: 

LOWER  LEFF  OR  VMAX,  UPPER  1-EFF  OR  VMAX,  STEP, 

LOUVER  [vGs],  UPPER  [vGs],  (sTEp) 

After  entering  these  last  values  plots  of  the  type  shown  in  Figs.  fS.l.  and  5.2.  will  be 
displayed  on  the  screen. 
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PROGRAM  2 


leO^DIM  Y0<3>,YK24>,A4<4),B4<4),P9<24> 

III  PRim  ^  2  FOR  IDis=F<y.1AX)‘ 

140  PRINT  "TYPEiMTO)HSUB«TOX<U))M<U)>LEFF<U>iUSB" 

150  INPUT  T$,U8,N0,T0,M,L0,U5 
160  T=1 

170  IF  T$="N'‘  THEN  190 
180  T=“l 

190  LO=L01:1.0E“4 
200  T0=T01:1.0E“4 
^10  l4=Wtl.0E-4 

220  PRINT  “U0,UCRIT|UEXP,UTRA,UMAX<11/S)" 

’40  PfilHT  "l6wER^LEFF*0R  UMAX» UPPER  LEFF  OR  UMAX, STEP" 

250  INPUT  m,W2,W3 
260  IF  T6=l  THEN  290 
270  Wl=Ult.l.0E-4 
280  W2=U2t 1 . 0E-4 
285  W3=W3tl.0E-4 
290  U=Utl60 

300  N=INT<<W2-U1)/M3)+1 

310  DIM  X(M),Y(10,H>  ^ 

320  PRINT  "LOWER  /UGS/, UPPER  /MGS/, /STEP/ 

330  INPUT  W4,W5,W6 

340  Hi=1.45E+10 

350  E0=3.45E-13 

360  E2=1.04E“12 

370  V0=0.02G 

380  Q=l.eE~19 

390  F2=Tt2*U0*LOG<M0/Hl) 

400  C0=EO/T0 

410  Q0=SQR<21:Q*E2tN0TF2CT> 


I  1\UUIW\|)  L. 


V ouli  L  I  nu^uy 


I 

tn 

I 


420 

430 

440 

450 

460 

470 

430 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 


GO=Q0/<C0*SQR<F2lcT)) 

FO-U3“F2-T1:G0*3QR(T5»;F2) 

Gmi0/C0 

U7=U3+G0*<SQR(T*(U5-»-F2)  )-SQR<  W2)  ) 

11  =  1 

FOR  I=W4  TO  W5  STEP  W6 
06=TTI 

U1=T*<U6-F0+U5> 

U2=T1:(05+F2> 

Al=4/3tG0 

B0=-2tUl 

D0=-02t2-4/3-tG0TSQR<02t3)+2TUUU2 
U3=O6-F0-F2+T.tG0t2/2.l:a-SQR<l+Tt.4/G01-2t.<U6~F0>>) 
U=U2TL0G  <U  I  .tE2/C0/  <  TT-  <U6“U7--U3TU3) ) ) 

U=U0tEXP<:U) 

IF  U<=U0  THEN  590 

U=U0 

Jl  =  l 

FOR  J=Wl  TO  W2  STEP  W3 

IF  Te=2  THEN  660 

L0=JT1.0E-4 

IF  I1>1  THEN  690 

X(Jl)=L0t 10000 

GO  TO  690 

U=Jtl0OO000 

IF  I1>1  THEM  690 

X(Jl>=Utl.0E-6 

U9=MTL0/U 

C1='’21:G0T.M9 

B1=B0-2TO9 

D1=D0+2TO9*U1 

A=~B1 

B=A1*C1-4*D1 

C=-D14:<A112-4*Bl)“Clt2 

R=-<At2>/3+B 


^uuii  u  I  nueu ; 


770  S=2*(A/3)t3-A1:B/3+C 

7S0  R3=R*R*R 

790  S2=StS 

800  P=S2/4+R3/27 

SIO  IF  P=>0  THEN  900 

820  P0=rtBS<P> 

830  RO=SQR(S2/4+P0) 

840  P2=-SQR(Pa> 

850  F=ATH(2TP2/3) 

860  R0=EXP(LOb<R0>/3) 

870  Y0(n=2TR0tCOS(F/3)-A/3 

880  N=3 

890  GO  TO  990 

900  H=1 

910  P2=SQR(P> 

920  P3=-S/2+P2 
930  P4=-S/2-P2 
940  P5=SGN<P3> 

950  Pb=SGN(P4> 

960  P3=P5TEXP<L0G(hBS(P3>>/3> 
970  P4=P6tEXP(L0G(ABS<P4))/3> 
980  Y0<l>=P3+P4-A/3 
990  M1=0 
1000  K>1 

1010  A3=Alt2/4-Bl+Y0<K> 

1020  B3=Y0(K>T2/4-Dl 
1030  A3=SQR<A3> 

1040  B3=SQR<B3> 

1050  FOR  f1=l  TO  4 
1060  Sl=IHT<(M-l)/2) 

1070  A4(ri)=Al/2+-11<f1-l)*A3 
1080  B4(M>=Y0<K>/2+-ltSl*B3 
1090  D4=A4(M>TA4<ri>/4-B4<M) 
1100  IF  D4<0  THEN  1170 
1110  t11=Ml  +  l 


rKuiiKMi'')  c  (continued} 


1120  Yl<Ml>«-A4(M)/2-»-SQR<D4) 

1130  P9<Ml)«Yl(Ml)t4+A11tYl<m)t3+Bl*Yl(ril>t2+Cl*Yl(Ml)+Dl 

1140  M1=M1+1 

1150  Y 1  < M 1  >  =- A4  ( M ) .''2-SQR  C  04  > 

1160  P9(M1  >=Yl(Ml)t4+All:Yl(rH>t3-fBl*Yl<Ml)t2+Cl*Yl<Ml>+Dl 

1170  MEXT  M 
1180  M2=0 

1190  FOR  K=1  TO  Ml 

1200  IF  ABS<P9<K>>>1.0E-5  THEN  1270 

1210  IF  Y1<K><0  THEN  1270 

1220  M2=M2+1 

1230  IF  M2>1  THEN  1250 

1240  Y2=Y1<K> 

1250  IF  Y2<Y1<K>  THEN  1270 
1260  Y2=Y1<K) 

1270  NEXT  K 
1280  04=TTY2TY2-li5-F2 
1290  IF  T6=2  THEN  1320 
1300  Y<I1, J1>=TTU4 
1310  GO  TO  1330 

1320  Y<I1, Jl>=W/L0*UtC01:<<U6-F0“F2-U4^2)*M4-2/3|:G01:(Y2t3-SQR(U2>t3 
1330  Jl=Jl+l 
1340  NEXT  J 
1350  11=11+1 
1360  NEXT  I 
1370  11=11-1 
1380  J1=J1-1 
1  <:•■?=*=; 

1400  S0=LGT<Y<I1,J1)) 

1410  S1=INT<S0) 

1420  IF  S0-S1=>LGT(7.5>  THEN  1510 
1430  IF  S0-S1>LGT<5)  THEN  1490 
1440  IF  £0-Sl>LGT<2.5)  THEN  1470 
1450  £2=2.5 
1460  GO  TO  1530 


PROGRAM  2  (continued) 


1470  S2=5 
1480  GO  TO  1530 
1490  S2=7.5 
1500  GO  TO  1530 
1510  S2=10 
1520  S3=10 
1530  SO=32t:l0tSl 
1540  DIM  Y7<J1) 

1550  PAGE 
15b0  GOSUB  1650 
1570  FOR  1=1  TO  II 
1580  FOR  J=1  TO  J1 
1590  Y7(J>=Y(I,J> 

1608  NEXT  J 

1610  MOUE  X<l>,Y7<:i> 

1620  DRAW  X,Y7 
1630  NEXT  I 
1640  GO  TO  2040 

1650  REM  ttt  DRAW  AND  LABEL  AXES  ttt 
1660  UIEWPORT  25,120,20,90 
1670  A=XC1> 

1680  B=X(J1> 

1690  C=0 
1700  D=S0 

1710  TS=S2/4tl0t<Sl“l) 

1720  T9=0.02*(B-A> 

1730  WINDOW  A-T9,B,C-T8,D 

1740  T=T8 

1750  MOUE  A,0 

1760  FOR  1=1  TO  J1 

1770  DRAW  A+n“l>#(B-A)/<Jl-l),0 

1780  PRINT  ••bJJ";X<I) 

1790  MOUE  A+<I-l)*<B-A>/<Jl-n,C-T 
1800  DRAW  A+(I-l>T<B>A)/<Jl-l),C+T 
1810  MOUE  A+(I-1>*<B-A>/<J1-1),C 


1820  NEXT  I 

1830  T-T9 

1840  NOME  0,0 

1850  FOR  1=1  TO  S3+1 

f?-?  A,  (  I-l  >.t(D-C)/S3 

IbfO  PRINT  ''HddbO'' ; S2/S3l!f I-l ) 

18b0  NOME  A-T,  C  I-l  >T<D-C>/'S3 
Ib^O  DRAW  A+T,  i  I-l  >.t:i:D-C)/'S3 
1 9uo  MOUE  A , ( I - 1 ) t < D-C > /S3 
1810  NEXT  I 
1920  HOWE  A,C+2''3*S0 
1930  IF  T6=2  THEN  1960 

l-a  rn'r^  i'pyyuUyfUH^’UJf'HJSUJflbJT* 

1980  MOUE  A+l/3t»  B-A'>  ,C 
1990  IF  T6=2  THEN  20^ 

2040  end 


Appendix  3 

Example  circuit  SPICE2.G  Input 


SECTION  or  CMOS  LSI  CIRCUIT  FOR  SPICE2,G 


■OPTIONS  ABST0L=1U  ACCT  VNTOL=100N  NODE 
*  2  IN  NAND  DEFINITION 
■SUBCKT  NAND  12  3  4  5 

Ml  3  1  6  5N1W=40U  L=8U  AD=800P  AS=LOOP  NRD=2  NRS=1 

M2  6  2  5  5  N1  W=40U  L=8U  AD=280P  AS=2080P  NRD=1  NRS=6 

M3  3  1  4  4P1W=20U  L=8U  AD=335P  AS=1164P  NRD=1  NRS=3 

M4  3  2  4  4Pliy=20U  L=8U  AD=336P  AS=lie4P  NRO=2  NRS=8 

■ENDS  NAND 


•  INVTRTEIRS  210  AND211 

•  SUBCKT  INV210  12  3  4 

M5  2  1  4  4N1W=24U  L=8U  AD=828P  AS=1024P  NRD=40  NRS=17 

M6  2  1  3  3P1>Y=24U  L=8U  AD=1432P  AS=4040P  NRD=5  NRS=18 

•ENDS  INV210 

•  SUBCKT  INV211  12  3  4 

M7  2  1  4  4N1W=8U  L=8U  AD=200P  AS=512QP  NRD=7  NRS=65 

M8  2  1  3  3P1W=8U  L=8U  AD=160P  AS=404CP  NRD=3  NRS=B 

•ENDS  INV211 

•  3  IN  NOR  GATE 

•  SUBCKT  NOR  1  2  3  4  5  6 

MIO  4  1  6  6N1W=20U  L=BU  AD=1150P  AS=5120P  NRD=22  NRS=65 
+  PD=125U  PS=230U 

Mil  4  2  6  6N1W=20U  L=BU  AD=1150P  AS=5120P  NRD=22  NRS=65 
+  PD=125U  PS=230U 


M12  4  3  6  6N1W=20U  L=8U  AD=1150P  AS=5120P  NRD=75  NRS=65 

+  PD=125U  PS=230U 

M13  7  3  5  5P1W=80U  L=8U  AD=420P  AS=960P  NRD=-5  NRS=1 

M14  8  2  7  5P1W=60U  L=8U  AD=1020P  AS=420P  NRD=1  NRS=^5 

M15  4  1  e  5P1W=60U  L=8U  AD=960P  AS=1020P  NRD=1  NRS=1 

•  ENDS  NOR 

•  OUTPUT  BUFFER 

•  SUBCKT  OUTBUF  1  2  3  4  5  6  7  8 

M20  4  8  6  6  N1  W=384U  L=8U  AD=4146P  AS=3043P  NRD=5  NRS=-5 

+  PD=445U  PS=15101JOFF 

M21  4  8  5  5P1  W=824U  L=8U  AD=11560PAS=93o0PNRD=-5NRS=^5 
+  PD=1190U  PS  =  :.65GU 

M22  8  7  6  6  N1  W=128UL=8U  AD=1900P  A3=3048P  NRD=2  NRS=5 
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+  PD=210U  PS=1510U 

M23  7  8  6  6N2W=8U  L=40U  AD=i800P  AS=6048P  NPD=10  NRS=15 
+  PD=150U  PS=1510U  OFF 

M24  7  8  5  5P2W=8U  L=40U  AD=1000P  AS=9360P  N^D=6  NRS=12 
+  PD=150U  PS=1650U 

M25  8  7  5  5  PI  W=12BU  L=BU  AD=1740P  AS=9360P  NRD=4  NRS=4 
+  PD=205U  PS=1650U  OFF 

M26  1  2  7  5P1W=24U  L=0U  AD=1400P  AS=300P  NRD=3  NRS=3 
+  PD=175U  PS=110U  OFF 

M27  1  3  7  6N1W=24U  L=8U  AD=600P  AS=540P  NRD=14  NRS=24 

+  PD=245U  PS=200UOFF 

CFF  7  0  B4F 

.ENDS  OUTBUF 

*  CIRCUIT  DESCRIPTION 

CALFIO  1  0  80F 

R1  1  2  2.2K 

C2  2  0  3.3P 

R2  23  235 

CALFI2  3  0  384F 

R3  3  4  1,27K 

CALFI3  4  0  116F 

XinAND  4  5  6  20  0  NAND 

VFl  5  0  10 

C6  6  0  200F 

XNOR  8  6  7  920  0  NOR 

XINVl  9  10  20  0INV211 

CIO  10  0  67.5F 

VR36  7  0 

XINVO  11  8  20  0INV210 
VCSO  11  0  10 
C8  80  65F 

XOUT  8  10  9  13  20  0  15  14  OUTBUF 
CL  13  0  20P 

VIN  1  0  PULSE(0  10  ION  ION  ION  20 ON  40 ON) 

VDD  20  0  10 

.NODESETV(13)=10V(14)=OV(15)  =  10 
•TRAN  IN  lOON 

•PLOT  TRAN  V(l)  V(2)  V(3)  V(4)  V(6)  V(9)  V(10)  V(l3)  (0,10) 

•MODEL  PI  PMOS(TOX=0.125U  NSUB=3.2E15  VTO=-0.65  XJ  =  1.5U  LD=1. 125U 
+  JS=7,75E-5  PB=0,86  UO=191  UCRIT=4.14E4  UENR=0.18  UTR4=0.25 
+  CGB=3E-10  CGD=3.26E-10  CG3=3.26E-10  CJSlf=  IF ’'A'LAX=5E4  NEFF=3 
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+  RSH=14LEVEL=2) 

•  MODEL  N1  NM03(T0X=Q.125U  NSUB=10E15  VTO=0.924  XJ=1.9U  LD=1.5U 
+  JS=L24E-4  PB=Q.89  UO=622  UCR1T=3.97E4  UEXP=0,08  UTRA=0.25 
+  CGB=3E-10  CGD=4.35E-10  CGS=4.35E-10  CJSW=2F  VMAX=5E4  NEFr=3 
+  RSH=3  LEVEL=2) 

•MODEL  P2  PMOS(T0X=0-  125U  NSUB=3-2E15  VTO=-0^65  XJ=L5U  LD=1-  125U 
+  JS=7-75E-5  PB=0^86  UO=222  UCRIT=5^63E4  UEXP=0-15  UTRA=0^25 
+  CGB=3E-10  CGD=3-26E-10  CGS=3-26E-10  CJSW=1F  VMAX=5E4  NEFF=3 
+  RSH=14  LEVEL=2) 

•MODEL  N2  NMOS(TOX=0^  125U  NSUB=10E15  VTO=0^924  XJ=1-9U  LD=1-5U 
+  JS=1.24E-4  P3=0-89  U0=814  UCR1T=6^2E4  UEXP=0.13  UTRA=0-25 
+  CGB=3E-11  CGD=4-35E-10  CGS=4,35E-10  CJ3W=2F  VMAX=5E4  NEFF=3 
+  RSK=3  LEVEL=2) 

•END 
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